Three dimensional theory for light matter interaction 
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We present a full quantum mechanical three dimensional theory describing an electromagnetic field 
interacting with an ensemble of identical atoms. The theory is constructed such that it describes 
recent experiments on light-matter quantum interfaces, where the quantum fluctuations of light are 
mapped onto the atoms and back onto light. We show that the interaction of the light with the 
atoms may be separated into a mean effect of the ensemble and a deviation from the mean. The 
mean effect of the interaction effectively give rise to an index of refraction of the gas. We formally 
change to a dressed state picture, where the light modes are solutions to the diffraction problem, 
and develop a perturbative expansion in the fluctuations. The fluctuations are due to quantum 
fluctuations as well as the random positions of the atoms. In this perturbative expansion we show 
how the quantum fluctuations are mapped between atoms and light while the random positioning 
of the atoms give rise to decay due to spontaneous emission. Furthermore we identify limits, where 
the full three dimensional theory reduce to the one dimensional theory typically used to describe 
the interaction. 



I. INTRODUCTION 

For several applications in quantum information sci- 
ence, such as long distance quantum communication 
it is essential to create an interface linking the photonic 
states used for transmitting quantum information to a 
material state suitable for storing and processing the in- 
formation. The generation of the required strong coher- 
ent coupling of light to a single emitter has proven diffi- 
cult to achieve inpractise, although substantial progress 
has been made [J H, 0, H, @] • In recent years optically 
dense atomic ensembles has e merg ed as a promi sing al- 
ternative 0, 1 1, M El M El El El El El M MM 

l2l"l . |22| . In this approach one can for instance use clas- 
sical laser pulses to engineer a suitable interaction such 
that an incoming light field is reversibly stored into the 
coherence between, e.g., two stable ground states in the 
atoms @. 

Some experiments on atomic ensembles uses atoms 
that are enclosed inside a cavity to enhance the cou- 
pling [l8| . In this situation the cavity defines a unique 
mode of the light field and the theoretical description con- 
sists of describing a single optical mode coupled to the 
atomic ensembles. Most experiments are, however, per- 
formed with atoms in free space not enclosed in a cavity, 
and in this situation the theoretical description is more 
complicated. Typically this situation is described in a 
one dimensional approximation, where one only consid- 
ers a single transverse mode and solves a one dimensional 
propagation equation for this mode @, EH Ell • 

In this paper we explore the range of validity of the 
one-dimensional approximation by making a full three 
dimensional description of the interaction between light 
and an atomic ensemble. Our calculations directly ap- 
ply to an experimental situations similar to the ones de- 
scribed in Refs. [1, H E1E3, where the light is detuned 
far from the atomic transition, but we expect the gen- 
eral features of our results to be valid for a much broader 



class of problems. 

Some justification for the one-dimensional description 
may be found in the literature on superflouressence, e.g. 
Refs. [H [H H^|- In this context it was found that 
the one-dimensional description is valid provided that 
the Fresnel number is of order unity T = A/ XL rj 1, 
where A is the transverse beam area, A is the wavelength 
of the light, and L is the length of the ensemble. Based 
on this work it has been argued that it is also neces- 
sary to have a Fresnel number of order unity in order 
for the one-dimensional approximation to be applicable 
to the quantum interfaces between light and atomic en- 
sembles @, EH EH- It is, however, essential to realize 
that the physical situations are very different in the two 
cases. The work on superflouressence typically concerns 
the temporal distribution of the output light measured 
by impinging the outgoing light on a photodetector. Be- 
cause the photodetector just measures the incoming flux 
/, this is essentially a multi-mode measurement 
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where the the sum is over all modes m hitting the detec- 
tor, and each of these modes are described by the photon 



creation (annihilation) operators a. 



t ). In particular 



the sum here includes all transverse modes. This is in 
contrast to the quantum interface work, where one is in- 
terested in the out goin g state of a single light mode, e.g., 
in Refs. d H [Tol llll| the measurement is essentially a 
homodyne measurement of a single mode, defined by the 
field of the strong classical laser. In other experiments the 
outgoing light is sent through a single mode optical fiber, 
which filters out everything except a single transverse 
mode. Furthermore the superflouressence work applies 
to a nonperterbative situation with a large optical gain, 
whereas the quantum interfaces typically operates in the 
few excitation regime. The previous analysis is thus not 
applicable to the present situation and it is therefore not 
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to be expected that the condition T ~ 1 is the right 
condition for the validity of the one-dimensional approx- 
imation. In fact, the experiments in Refs. [H, H, E3] 
are performed with T ~ 10 4 , and still give very good 
agreement with the one-dimensional description. Here 
we make a full three dimensional description of the ex- 
periments in Refs. [8, §, and we find that it reduces 
to the one-dimensional description in the paraxial ap- 
proximation provided that T 3> 1. 

In a related work a three dimensional description was 
also presented in Ref. J26[. Whereas our procedure as- 
sumes non-moving atoms, i.e., cold atoms, that work con- 
sidered the opposite limit, where the motion of the atoms 
wash out any spatial structure of the atomic spin state. 
Unlike the situation in Ref. [26[ , where the motion of the 
atoms always lead to certain inefficiencies, the fact that 
we consider stationary atoms, allows us to identify cer- 
tain limits, where we exactly reproduce the simple result 
of the one dimensional theory as discussed in Sec. IVI Bl 

Our theory is developed as a perturbative expansion of 
the interaction between light and the atomic ensembles. 
It is, however, essential to be very careful about the way 
this perturbative expansion is performed. Below we shall 
present results up to second order in the interaction be- 
tween the light and the atoms. We shall use an effective 
Hamiltonian, where the excited atomic state has been 
eliminated, i.e., a Hamiltonian of the form 

H ~ ^2Y1 5k,k'"k(ri)ttk/ (r;)a k ,d k , (L2) 

k,k' i 

where <?k,k' is a coupling constant for the two modes k, 
and k' described by photon creation (annihilation) op- 
erators d k (dk) with mode functions Uk, and is the 
position of the ith atom. If we take the mode functions 
to be simple plane waves with an input field in a certain 
mode ko and calculate the intensity in a certain direction 
described by ki, we find the intensity 
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(1.3) 



where Ak = ki — ko. The standard way to proceed 
from here is to say that the exponential varies rapidly 
when i ^ j and therefore neglect all terms except i = j 
so that one is left with something proportional to the 
number of atoms Na, which is known as spontaneous 
emission. For the problem we are interested in here, we 
are, however, mainly concerned with the properties of 
the light in the forward direction, where Ak sa 0. In this 
case it seems unjustified to neglect the cross terms which 
give rise to collective scattering scaling as iV 2 . Since N is 
typically a very big number, the presence of such large N 2 
contributions may limit the applicability of perturbation 
theory. 

In order to avoid the problems associated with this 
collective scattering, we use a different basis for our per- 
turbative expansion: instead of starting from the eigen- 
modes of the propagation equation in vacuum, we use 



the solutions to the classical diffraction problem in the 
presence of the medium, i.e., we take into account that 
the atoms give rise to an index of refraction of the gas, 
which changes the propagation of the light. Specifically, 
we write the Hamiltonian as 



H = (H) a 



■SH, 



(1.4) 



where (7i) atoms is the quantum mechanical expectation 
value of the Hamiltonian with respect to the atomic 
spin state averaged over the random positions of the 
atoms. This averaged Hamiltonian gives rise a contin- 
uous quadratic Hamiltonian in the light field operators 
similar to a Hamiltonian describing the interaction with 
a dielectric medium. When we formally change to the 
interaction picture with respect to this averaged Hamil- 
tonian, we obtain a new set of basis modes. Doing per- 
turbation theory on these modes, the only effect on the 
light comes from the quantum mechanical fluctuations 
and the fluctuations caused by the random position of the 
atoms. These fluctuations are described by the Hamilto- 
nian SH = H — (H) atoms- When we average the first order 
term in the perturbative expansion with respect to the 
position of the atoms the resultant expression describe 
that the quantum fluctuations of the atoms are mapped 
onto the light in analogy with the results derived in a 
one-dimensional theory in Ref. Q ■ 

If we go to second order in the interaction, our expres- 
sion will give terms quadratic in SH. In order to take 
the spatial average of such terms we need to know the 
density-density correlation function of the atoms. In- 
serting the density-density correlation function for an 
ideal gas we no longer find the collective scattering terms 
described above, i.e., the collective scattering is essen- 
tially the classical diffraction of the light, which is ex- 
plicitly taken into account by our average Hamiltonian, 
and therefore it does not appear in our perturbation the- 
ory. The spatial average of the second order term does, 
however, produce a new term associated with the point 
particle nature of the atoms and their random positions. 
This term is equivalent to the results obtained by just 
keeping the i = j terms in Eq. (|I.3p . and represents the 
effect of spontaneous emission. 

Unlike most approaches to the interaction between 
atoms and light, which derive coupled equations for the 
atomic states and the electric field, our approach consid- 
ers the electric displacement field D instead of the electric 
field. The reason we chose to use the displacement field 
is that it is convenient to work with a purely transverse 
field, which is the case for the displacement field due to 
the macroscopic Maxwell equation V ■ D = 0, whereas 
this is not necessarily the case for the electric field in a 
medium. Formally the two approaches are equivalent and 
may be related through a unitary transformation [27| . 

The full theory is quite involved. Readers who are 
mainly interested in the consequences of our theory for 
experimental implementations are therefore advised to 
skip to Sec. IVI[ where we discuss such consequences. The 
sections prior to this mainly focus on building the theo- 
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retical frame using a first-principles strategy. The paper 
is organized as follows: In Sec. [II] we give the details of 
the model used to describe the interaction. In Sec. IIIIIwc 
derive a set of equations of motion describing the system 
of atoms and light, using Heisenberg's equation of mo- 
tion. The wave equation describing the light is expressed 
in a form that ideally suits a perturbative treatment. In 
Sec. HVI we express the general solution to the wave equa- 
tion in terms of Green's functions and derive the pertur- 
bative expansion of the solution to the wave equations 
as well as the equation describing the atoms. This is 
represented in terms of Feynman diagrams. In addition 
we develop the appropriate theoretical tools to describe 
point particle effects such as density correlations, and de- 
rive a formal expression for the Green's function. In Sec. 
IVl we present our results where we discuss higher order 
effects such as spin decay and light scattering. We de- 
fine operators that describe photon-measurements, and 
demonstrate how these are calculated in the theory In 
Sec. I VII we discuss various limits where the general three 
dimensional theory reduce to the usually employed one 
dimensional model We also describe how a detailed 
understanding of the spatial modes can be used to achieve 
storage and retrieval of information in several transverse 
modes of light and atoms simultaneously. In Sec. I VIII we 
conclude the work, and in the appendices we give several 
details omitted from the main text. 



I F",m > 



I F',m+1 > I F',m > I F',m-1 > 




I F,m+1 > I F,m > I F,m-1> 

Figure 1: Example of an atomic level structure. The atoms 
have a single ground level with spin F and one or more ex- 
ited levels. The fields have a large detuning A so that the 
exited states may be adiabatically eliminated and we obtain 
an effective ground state Hamiltonian Eq. (|II.3|I . 



the polarization of the media [27|, |26 



Atoms 



H 



int = - ]T — D(rj,t) • P(rj, f). 

A 6 



(11.1) 



II. MODEL 

The model we consider describes the interaction be- 
tween an ensemble of atoms and an incoming light field. 
The atomic ensemble is considered to be an ideal gas of 
identical atoms. The atoms are described as non- moving 
randomly distributed point particles and the interac- 
tion with the light field is described within the dipole- 
approximation. Each atom is assumed to have a ground 
level of total spin F. In addition we assume that the 
atoms have no other stable ground states to which they 
can decay. See Fig. [I] We shall assume that the elec- 
tric fields are sufficiently far-detuned that we may adia- 
batically eliminate the exited states, and work with an 
effective Hamiltonian involving only the ground states. 
In the following we first discuss the interaction between 
light and a single atom, and then move on to discuss the 
interaction with an ensemble of atoms. 



A. Interaction with single atoms 

The aim of this work is to describe the interaction be- 
tween an electromagnetic field and an ensemble of iden- 
tical atoms. The problem is therefore both to deal with 
the microscopic behaviour of a single atom, and also the 
collective effect of many atoms. We choose here to work 
in the so called length gauge, where the basic interaction 
is given as the product of the displaced electric field and 



Our gauge choice ensures V • D(r, t) = 0. We will assume 
that the fields have a large detuning and do not saturate 
the atomic transition, so that the exited levels may be 
adiabatically eliminated. This procedure is described in 
Appendix [X] The polarization of the atomic ensemble 
then depends linearly on the displaced electric field, that 
is P(r, t) = I/[J]D(r, t). We introduce here the argument 
J to indicate that the interaction matrix V[3] depends 
on the spin of the atoms. Next we write the displaced 
electric field as a sum of a positively oscillating part and 
a negatively oscillating part, 



D(r,i) = D (+) (r,i) +D ( ~ ) (r,t). 



(II.2) 



In Appendix fX] we show that the effective interaction 
Hamiltonian, assuming such linear dependence of the po- 
larization on the displaced electric field, reads 



Atoms 



D 



(+) 



D 



(-) 



V[3 3 ] D 



(+) 

3 



(II.3) 



where we have also employed the rotating wave approx- 
imation. Here the superscript t denotes matrix transpo- 
sition. 

Since the Hamiltonian must be rotationally invariant it 
can only contain irreducible tensors of at most rank two. 
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In the vector representation the interaction may thus in 
general be written as 

?[3j] = p(co X} - ic x J j x +c 2 x Jj) ■ (Jj x ) . (II.4) 

The meaning of the notation is that when inserted into 
the Hamiltonian the result of, e.g., the last term of the 
right hand side of Eq. pi.4[) is 

Atoms 

/3c 2 £ (D(-)(r„t)xJ,)-(J, xD(+)(r„<)). (II.5) 

3 



Note that we have here chosen a description which has 
a simple analytical representation, but this means the c 2 
term is not a pure rank two irreducible tensor, but consist 
of a combination of tensors of rank zero, one and two. In 
matrix form the interaction may be written: 



V[3] = (3 



(co - c 2 )J 2 + c 2 J 2 ic x J z + c 2 J y J x -ici_Jy + c 2 J z J x 
-ic 1 J z + c 2 J x J y (c — c 2 ) J 2 + c 2 Jy ic x J x + c 2 J z J y ^ 
ici J y + c 2 J x J z -ici J x + c 2 J y J z (co - c 2 ) J 2 + c 2 J| 



(II.6) 



In general the atoms may have several exited levels as 
shown in Fig. [T] The effect of several exited levels can 
be included in the coefficients cq, c\ and c 2 that will then 
depend on the detuning. For atoms with F = ^ or for an 
alkali atom, where the fields are detuned by more than 
the hyperfine structure of the exited state, the c 2 term 
disappears p9l | and the interaction matrix is given by 



V[Jj] = /3 co 3^ -id J 



(II.7) 



Here Co and c\ are constants which depend on the atomic 
structure as well as the detuning. The coupling constant 
[3 in Eq. pi.7|) is given by 



(3 



7T7 



(II. 



where 7 is the linewidth of the exited level, A the de- 
tuning of the laser field with respect to the atomic tran- 
sition, and k L is the wave vector. With this choice of 
(3 the coefficients Co, c% and c 2 will be of order unity or 
less. Throughout this paper we shall only consider the 
simple interaction in (|II.7|) . A discussion of the effect of 
the c 2 term is given in Refs. [U H3] in a one dimensional 
description. 

We will consider a perturbative regime, where the 
product of the atomic density p and (3 is small f3p <C 1, 
and make a perturbative expansion in (3. Note, however, 
that this condition does not imply that the total effect 
of the interaction is small. On the contrary, we are most 
interested in situations, where the integrated effect of the 
interaction significantly alters the light beam as it passes 
through the sample. To take into account these collec- 
tive effects we explicitly include, e.g., the diffraction of 
the light caused by the propagation through a medium. 
To describe these effects we discuss in the following sec- 
tion how to quantize the field in a medium. 



B. Mode expansion 

To quantize the electromagnetic fields we could: i) im- 
pose the canonical commutation relations on the vector 
potential and displaced electric field. Or ii) expand the 
electromagnetic fields on an orthonormal set of spatial 
mode- functions {fk} conveniently chosen to diagonal- 
ize the Hamiltonian (in vacuum this is the set of plane 
waves), and then quantizing the mode-amplitudes. Here 
we will use the latter. The Hamiltonian describing the 
electromagnetic field in a medium is given by (27j 



H = 



-D 2 (VxA) 2 



Mo 



H 



in 1 3 



(II.9) 



where 7ii n t is given in equation ((II . 3 j) . A careful analysis 
of how to quantize the electromagnetic field in a medium, 
is given in Ref. (301, and here we shall only go through 
the steps briefly. 

By introducing the spin field 



t) 



Atoms 



5{r 



l 3)i 



(11.10) 



the Hamiltonian may be put in an all-integral form. The 
main idea in our approach is to divide the full Hamilto- 
nian into a spatially averaged part, and a point particle 
part, describing the fluctuations from the average caused 
by the atoms being point particles. For now we only 
consider the spatially averaged part of the theory. We 
will use calligraphic font to denote that we have made 
a spatial average. We thus write the spatially averaged 
interaction from equation (|II.7p as 



V[J] =/3 P (r) c J 2 -ici J(r) x 



(11.11) 



5 



Here a bar denotes a single-atom operator, that is J(r) 
is the spin operator of a single atom at position r. Wc 
use the bar to distinguish between the spatially aver- 
aged single-atom spin operator, and the general spin 
field in equation pi.lO[) . The two may be related by 
(J(r, i)) s . a . = p(r)J(r,t), where (-) s . a . denotes spatial av- 
erage. The function p(r) denotes the average atomic den- 
sity, which in this model is a continuous scalar field. 

In the following we will define a mean Hamiltonian, 
where we have taken into account the quantum mechani- 
cal average of the spatially averaged interaction. We then 
write the Hamiltonian as a sum of the average Hamilto- 
nian and a point particle Hamiltonian 



H =Hq + Tip 



(11.12) 



where 



H =- / d*r{ 



- D(.M'D(-) + MDM) (V x A) 



Mo 



■}• 



n 



— J d 3 r D • (fhlJf D ( - } + m[J] D (+) 



(11.13) 



|,] '~ 2e 

M = i- v[J] 



and 



ift[J] = V[J|-V[J]. 



(11.14) 
(11.15) 



(11.16) 



Here we simply write J (without the hat) to denote that 
this is now a classical field describing the classical expec- 
tation of the spin of the atoms. In analogy with Ref. [3(| 
we introduce the mode functions {fk} defined by: 



The minus sign in Eq. pi,19a[) is conventional and stems 
from the relation between the displaced electric field and 
the canonical conjugate field given in terms of the vector 
potential. 

The reality condition on the displaced electric field 



(D(r,t))t =D(r,t) 



allows us to write 



D(r,i) 



k 



(p+(t)f k (r)+p k (t)f£(r)). (11.20) 



Using the results in E qs. (|II.17[) and (|H-18|) and the 
expansion in equation pi,19|) . the Hamiltonian attains 
the desired diagonal form 



n =\ J d»r{ 



D(1-V[J])D (VxA) 



Mo 



\Y, {PIMM*) + "&£(*)<&(*)}• (11.21) 



The mode functions {fk} are thus the spatial basis di- 
agonalizing the spatially averaged Hamiltonian, and as 
we shall see the proper basis describing the diffraction 
problem. 

The splitting in equation (|II.12j) allows us to consider 
the problem as comprised of two types of properties. The 
effect of single atoms, and the spatially averaged Hamil- 
tonian. The effect of the spatially averaged Hamiltonian 
is well understood in terms of the mode-functions defined 
in equation (|II.17[) . The point particle effect we will dis- 
cuss in greater detail when considering the equations of 
motion for the full system. Before deriving these equa- 
tions of motion we, however, briefly need to discuss the 
commutation relations describing the system. 



V x V x Mf k (r) =-|fk(r), 
V • f k (r) =0. 



(II.17a) 
(II.17b) 



We also define the appropriate inner product on the space 
spanned by these mode functions: 



(0(r)|V(r)) 



d 3 r0(r)* -Mip{r). 



(11.18) 



We will assume that the average interaction term V[J] 
does not evolve in time, and our appropriate mode- 
functions are therefore time independent vector fields. 
One can show that the functions fk span a complete or- 
thonormal basis for the space in which we work. To diag- 
onalize the Hamiltonian we expand the vector potential 
and the displaced electric field in these mode functions 

r>(r,t)=-J2V^P^K(r) (II.19a) 

k 

A(r, t)=J2 cVJ^ «k(*)(l ~ ^[J])/kW. (IL19b) 



C. Quantization and Commutation relations 

Above we expanded the fields in convenient spatial 
modes. The coordinates Pk(i) and qk(t) are canonically 
conjugate variables, and we can thus quantize our theory 
by imposing the commutation relations 



[qu(t),pk'(t)] = ih& 



kk' 



(11.22) 



It will however be convenient to have the commuta- 
tion relations for the fields which we may derive from 
the mode-amplitude commutation relations. It will also 
be convenient to separate the displaced electric field 
into a positively and a negatively oscillating part D = 
D(+) + D(-), where D<-> is in accordance with conven- 
tion chosen so that it only contains terms oscillating like 
e luJkt . Our choice of gauge is reflected in the transver- 
sality of the mode functions defined in Eq. (|II.17|) . We 
expect this transversality condition to be represented in 
the commutation relations as well. With the quantiza- 
tion procedure above one finds the following expression 
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for the negative frequency part of the relevant fields 



^4e***5(r) (II.23a) 



A(-)(r,<)=^ C ,/|^4 e ^(l-V[J])^(r). 



(II.23b) 



The positive frequency part may be found by Hermitian 
conjugation. The above result is found from equation 
([TT720]) along with the definit ions of creation and annihi- 
lation operators given by 



ft(<)=\/^ : {« k (i)+E c/ k k '«Lw} 

k' 



Pk(t) =i 



(II.24a) 
(II.24b) 

(11.25) 



A detailed discussion of this procedure is found in Ref. 

From these definitions and the commutation relations 
pi.22j) we obtain 



where the matrix £/ kk / is defined as 

Pick' - / d 3 r^f k (r)-f k ,(r). 



[a k (i), &£(*)] = <5 k k'- 
Going to the field operators we get 
[D (+) (r,i),A (+) (!■',<)] =0 

[D(+)(r,f),A<-V,*)] =p T (r,r>) 



where 



I T (r,r')=^f k (r)[^(r 



(11.26) 

(11.27) 
(11.28) 

(11.29) 



Here <5 T (r, r') is a generalized transverse delta function 
[30| . This may be seen by considering its action on some 
transverse vector field (V • ip(r,t) — 0). Since {fk} is 
a complete basis on the set of transverse fields, we may 
expand V( r > as 



^(r,t)=^C k (t)f k (r). 



(11.30) 



If we calculate the effect of the transverse delta-function 
on a transverse field we find 



dV <5 T (r,r') • ip{r',t) = 



d 3 r' J2 C'kWfk'(r) \M%(r') ■ f k (r') 



kk' 



^(7 k (t)f k ,(r)5 kk , =V(r,f), 



(11.31) 



where we have used the orthonormality condition of the 
basis-functions. 

We shall also need the equal-space commutation rela- 
tions 

[D<+>(r,f),D<->(r,f')]. 

A formal expression of this commutation relation can be 
found from Eq. (|II.23a|) to be 

[E>M (r,i),E)(->(r,0] = ^V(r,t,t'), (IL32) 



where 



(r,t,tO=^a; k f k (r)^(r)e- 



-2CJk (t — t') 



(11.33) 



In vacuum ^(r, t, t') is simple to evaluate, but for complex 
systems it is nontrivial to gain knowledge of the basis- 
functions {f k }- In Appendix IBl we calculate fj using the 
rotating-wave approximation and the local density ap- 
proximation, where we assume that p(r) varies slowly 
with respect to r. 



III. EQUATIONS OF MOTION 

In this section we derive the equations of motion for 
the system, and consider their general properties. In the 
previous section we discussed that the theory could be 
divided into an average part and a part representing the 
deviation from the average. To derive the equations of 
motion we will, however, work with the full Hamiltonian 
and later make the splitting into the average part and 
the deviations from it. The strategy we will use is to 
first derive the quantum mechanical Maxwell equations, 
and then to combine them into an effective wave equation 
for the field. 

We will now as an example derive one of the quantum 
mechanical Maxwell equations from Heisenberg's equa- 
tion of motion: 

=^-J d " r ' [(Vx A(r')) 2 ,D(r)] 

J d 3 r> {(V x V x A(r')) • [A(r'),D(r)] 
+ [A(r'),D(r)] ■ (V x V x A(r'))}. (HI.l) 



Here we have used the Hamiltonian given in Eq. pi.9p . 
and the boundary condition that the physical fields van- 
ish at infinity. To shorten the notation we have sup- 
pressed the explicit time dependence. The commutation 
relation may be found from pT27| and (|H2"8)) to be 



kk' 



[A(r'),D(r)] = -ihF(r,r') 



(III.2) 
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Since the field V x A is transverse by definition, this 
gives us the first quantum mechanical Maxwell equation. 

^D(r) =-V x B(r), (III.3) 

at /iq 

where 

B(r) =V x A(r). (III.4) 

Similarly we may derive the Maxwell equation V x E = 
— 9tB, where E = —dA/dt = D P. The remaining 
Maxwell equations V • B = and V • D = follow 
immediately from the definition of B in Eq. (|III.4[) and 
from the transversality of D. 

Because of the nature of the interaction part of the 
Hamiltonian, it is convenient to consider the two fre- 
quency components of the displaced electric field sepa- 
rately. The quantum mechanical Maxwell equations may 
be combined into a single wave equation 

/ d 2 \ 

(d~2 +c2y x V x )V { - } (r,t) = 

c 2 J d 3 r V x V x <F(r,r') • f [Jftt^ (r' ,t), (III.5) 

where the positive frequency part can be found by Her- 
mitian conjugation. Similarly we may derive equations 
for the spin of the atoms, and for the simple interactions 
given in Eq. f|II.7[> . one finds 

p(r,t) = ^J(r,t) x (D(-)(r,t) x D<+>(r,t)). 

(III.6) 

In the remainder of this article we will solve these coupled 
partial differential equations. 

The expression in Eq. (lIII.5j) is a second order differen- 
tial equation in time. The solution of this equation will in 
general not only depend on the initial value D(r, t = to), 
but also the time derivative <9 t D(r, i)| t=to . In deriving 
our interaction we have, however, already used the rotat- 
ing wave approximation, where we ignore the dynamics 
on a time scale similar to the inverse of the optical fre- 
quency. Similarly we shall here make a slowly-varying- 
envelope approximation and write the displaced electric 
field as 

D(r,<) =D(-)(r,t)e lWLt +D (+) (r,i)e"^ Lt , (III.7) 

where are slowly varying in time. If we ignore 

the second derivative of the slowly varying operators 
(<9 2 DW(r,i) w 0), then Eq. (|nT3|) reduces to a first- 
order differential equation in time. 

Since we are heading towards a perturbation theory in 
the point-particle part of the Hamiltonian pi,12[) . we will 
add and subtract the average part of the source term in 
Eq. pil.5p . That is we write 

V"[J] = V[J] - V[J] + V[J] = fh[J] + V[J]. (III.8) 

The idea in this separation is that now V[J] represents 
the average effect of the ensemble, which may have a big 



effect, whereas to [J] represents the fluctuations around 
this average. To take advantage of this we first consider 
the average term 

J dVVxVxFtr.r'j'VlJfDHlr',!). (IH.9) 

This term is continuous and we may use partial integra- 
tion twice. Using the expression for the general trans- 
verse delta-function one finds 

J d 3 r V x V x I T (r,r') • V[J]'D(->(r', i) 

= VxVxV[JfDH(r,(). (111.10) 

This term we will move to the left hand side of Eq. 
pil.5|) . and we are left with a diffusion equation involving 
only the fluctuations as a source term on the right hand 
side 

(2iu^-ul + c 2 VxVxjM') t) { - ] (r, t) 

= c 2 J d 3 r V x V x l T (r,r') • m[J]*D(-> (r' ,*). 

(111.11) 

If we put the right hand side of this equation to zero, 
i.e., ignore the fluctuations, this equation describes the 
propagation and diffraction of the field in a medium. For 
instance if we take the simplest case where the medium 
is isotropic so that the matrix V[J] is just a scalar, this 
equation describes the propagation through a medium 

with an index of refraction given by n = 1/ \J 1 — V[J], 
seeRef. (H- 

IV. GENERAL SOLUTION AND FEYNMAN 
DIAGRAMS 

In this section we discuss the solution of Eq. (|III.11[) 
in terms of its Green's function. Let us for convenience 
define the differential operator 

V = 2iuj^ - luI + c 2 V x VxM*(r). (IV.l) 
We then define the Green's function by 

vd^(r,t\v ,t ) = F(r,r )5(t-t ). (IV.2) 

The right hand side of this equation describes an identity 
functional on the inner product space we are working in. 
We want the Green's function to describe an evolution of 
the system forward in time. We therefore define a cut-off 
on the Green's function in time 

&(-\r,t\r ,t ) =0 for t<t . (IV.3) 

The general solution to Eq. pil.ll[) in terms of Green's 
functions is discussed in detail in Appendix [Cl and reads 
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D(")(r,f) =21^ \ d 3 r' M t (r')d^(r,t\r',t )-t>(-\v',t ) 

fr / / d 3 r'dt' ^t*(r')(5(-)(r,t|r',0 • / dV' V x V x l T (r',r") • m[3]*D<->(rV)- 0V.4) 



The upper limit is understood to be i - *" — lim e ^o [t ~\~ s] ■ 
Before continuing a few comments are in order. Here we 
have used the boundary conditions, that all fields van- 
ish at infinity, i.e., we imagine that at time t — we 
have generated an optical pulse inside the volume we are 
describing, which travels toward the atomic medium. Al- 
ternatively we could have described the incomming field 
by a boundary term. The positive frequency part may 
be found by Hermitian conjugation. 

Let us now consider the last term of Eq. (|IV.4j) . We no- 
tice that the involved fields are all continuous and differ- 
entiable with respect to the primed spatial coordinates. 
Using partial integration twice and introducing the prop- 
agator defined by 

p(-)(r,*|r',i') = V x V x .M'tOG^MIiV) 

(IV.5) 

the last term of Eq. PV.4P may be written as 

c 2 ff d 3 r'dt' fd 3 r" p(-)(r,f|r^><F(r',r") 

•7fi[J]*D(-)(r",i')- 
(IV.6) 

Due to the cross product in Eq. pV.5[) the propagator 
is transverse with respect to primed coordinates and the 
transverse delta function in (|IV.6p may be integrated out, 
giving 

c 2 d 3 r'dt' p(- ) (r,t\r',t')-fh[J] t f) { - ) (r , ,t'). (IV.7) 

The first term of the right hand side of equation (jIV.4l) 
we will denote as Dq (r, f) 

D^(r,t) = 

2ilu l J d 3 r' M*(r')GM(r,i|r',i ) . f)^(r',t ). 

(IV.8) 

If there were no deviation from the mean, i.e. m[J] = 0, 
the solution would simply be D^^(r,t) = '(r,t). 



Dq (r, t) thus denotes the solution to the diffraction 
problem, where the atomic medium is treated as a con- 
tinuous medium with a diffraction matrix M.. 

A. Perturbative expansion 

Below we shall develop a perturbative expansion in the 
deviation from the mean due to quantum fluctuations 
and from the fact that the medium is not continuous 
but consists of a large number of point particles. The 
starting point for the perturbative expansion will be the 
field equation 

D(-)(r,t)=D<- ) (r,t) 

+ c 2 J J d 3 r'dt' P^(r,t\r',t')-m[J] t 'b^(r',t'). 

(IV.9) 

In addition to this we shall also need the solution to the 
equations of motion for the spin (|III.6|) . which may be 
formally solved to give the spin equation 

J(r,t)=J(r,to) 

+ l l2L I dt > j( r ,i') x (D(-)(r,<) xDWfr.i')). 
J t v ' 

(IV.10) 

These are the equations we wish to treat using the Born 
approximation, where we make an expansion in the in- 
teraction parameter 0. ( In Eq. (|iV.9[) the interaction 
m[J] is proportional to the expansion parameter (3.) 

In terms of notation this expansion gets extremely 
cumbersome. It is therefore convenient to introduce 
Feynman diagrams to represent the various terms of the 
expansion. We will be dealing with two types of interac- 
tions: the one given in Eq. |E9]) which we will represent 
with a shaded circle, and the one given in Eq. (jIV.lOp 
which we will represent with a shaded triangle. The field 
equation, we diagrammatically represent as 




(IV. 11) 



and the spin equation is represented as 



(IV.12) 



The orientation of the diagram is such that time is going 
from left to right, and the evaluation at time t is marked 
by a dot. Spin propagation is represented by a line with 
an arrow pointing in the positive-time direction. A wig- 
gly line represents propagation of the displaced electric 
field. The arrow denotes whether the line represent the 
photon-generating part of the field, )(r, t), where the 
arrow points forward in time, or the photon-annihilating 
part of the field, where the arrow points backward in 
time. The full solution to the spin 3(t) is denoted with a 
double straight line, and the full solution to the displaced 
electric field is denoted with a double wiggly line. 

The field equation and the spin equation can be repre- 
sented as a perturbation series, and in the following we 
shall discuss the effect of the terms in this perturbation 
series. An important feature of our system is the random 
distribution of the atoms in the ensemble. The equations 
that we have derived so far apply to each realization of 
the atomic distribution {ri, r2, . . . , rjv} ■ However since 
we have no control of the position of the atoms we will 
have to make a spatial average of our equations, that is 
of the terms in the perturbation series. To do this we 
need to know the density correlations of the gas. 



B. Density correlations. 



We assume that we are dealing with an ideal gas, i.e., 
we assume that the distribution of the atoms is com- 
pletely random but has a distribution given by the pos- 
sible spatially varying density p(r), and we assume that 
there are no correlations between the positions of dif- 
ferent atoms. The correlation function for the density 



distribution p(r) — £V H r ~ r j) ^ s thus 
<p(r)p(r'))s.a. = (E^-^Wr'-r ; )),a. 

= ^(<5(r-r J )<5(r'-r i )) s ., 

+ J2S(r-r'){S(r-r j )) s . a ,. 

3 

= (p(r)) s . a .(p(r')) 8 . a .+^(r-r')(p(r)) 8 . a .. 

(IV. 13) 

Here (-^.a. denotes spatial averaging. In the last step 
we used that the distribution is independent for different 
atoms, and we ignored the small difference between N\ 
and Na{Na — 1), where Na is the number of atoms. We 
have also neglected the effect that two different atoms 
can not be found at the same point in space. While this 
may seem insignificant for a low density gas, we show in 
Appendix [D] that including this effect to all orders in the 
perturbation series gives the Lorentz-Lorenz correction 
to the index of refraction. 

Below we shall also use the correlation functions for 
the spin. Similar to the calculation above we find 

(J„(r)J m (r')) s . a . = p(r)p(r') J„(r) J m (r') 

+ p(r)<5(r-r / )J„(r)J m (r), 

(IV. 14) 

where the index n, m refer to the spatial components of 
the operators. To shorten notation we have written p(r) 
instead of (p(r)) s , a ,. As discussed previously the bar de- 
notes a single atom operator. We will preserve the quan- 
tum mechanical behavior of the operators by not taking 
the quantum mechanical mean. The first term on the 
right hand side of Eq. pV.14[) arises from the contribu- 
tion from different atoms (signified by the prime on the 
second spin operator). In the second term on the other 
hand the two operators refer to the same atom, and the 
operator product should be evaluated for a single atom. 
For example for a spin-^ system, we have the follow- 
ing relation between products of spin operators on single 
atoms 



J«(r)J m (r) 



(IV.15) 
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The generalization to even higher-order density correla- 
tions is straight-forward. 

These considerations become important when we cal- 
culate the spatial average of the second-order terms of 
the perturbation series. Let us as an example consider 
the second order term of the spin equation representing 
a photon first interacting with one atom and then later 
with the atom in consideration. 



s.a. 




(IV. 16) 




When taking spatial average this term generates two 
terms in the perturbative expansion as indicated with 
the arrow in Eq. pV.16p . The first term involving the 
spin of two different atoms we will refer to as a coher- 
ent interaction, which we will discuss later. The second 
term involving the delta function corresponds to the in- 
coherent interaction (for reasons which will become clear 
below). We include this situation in the diagrammatic 
notation by introducing a hatched star and a loop sig- 
nifying the infinitely short propagation stemming from 
the delta-function term of the correlation function Eq. 
figjl) , i.e. 



d?rP ± (r,t\r',t')-i/)(r',t')5(r-r') = 



P ± (r,t\r,t')-if,(r,t'). 

(IV. 17) 

The loop is placed on the top of the star when it comes 
from the positively oscillating propagator p(~\ and in 
the bottom of the star when we refer to the negatively 
oscillating propagator p( + \ A star scales with the ex- 
pansion coefficient (3 squared since it involves two interac- 
tions. In the next section we will calculate the infinitely 
short propagator appearing in these expressions in the 
local density approximation. 



C. Green's function and propagator 

In this section we first derive a formal expression for 
the Green's function. Within our inner product space 
the Green's function is defined by (|TVT|) and (|TV2|) . Ex- 
panding our Green's function in the basis f^(r) we find 
the representation 

<3<-> (r, £|r', = f k (r)fk(r')^- } (t, f )• (W.18) 



We have here expanded on the complex conjugated set 
fk(r) to match the expansion of the displaced electric 
field in Eq. (|II.23aj) . The transverse delta- function has 
the representation 



I T (r,r') = ^f^(r)f k (r') 



(IV.19) 



where we are now working in the inner-product space 
with inner product defined in Eq. (|II.18j) . The scalar 
function g^\t,t') is defined by 

(2iu^ - ul + u£)gt\t,t') = S(t - t% (IV.20) 

along with the condition that the function g k (i, 1') vanish 
for t < t' . We will consider the following form of the 
scalar function, where we explicitly write this cut-off in 
terms of a step function 



■9k 



{ -\t,t') = Ce i ^ t - t '^Q{t-t l ). 



(IV.21) 



The coefficients 7k and C is found by inserting this result 
into equation (|IV.20j) . 



7k : 
C 



_ ^k~^ 
2w L 
—i 
"2^7' 



The Green's function is thus given by 
G(-)(r,i|r',0--*E f k( r )fk(O- 



2w L 



(IV.22a) 
(IV.22b) 

-e(i-t'). 

(IV.23) 



Next we will look at the infinitely short propagator in 
Eq. pV.17p . Using the Green's function given in equa- 
tion (IIV.23P along with definition pi.!7a|) and pV3|) the 
propagator may be written as 



#(-)(r,t|r,i') 



2lu,c 



^^(r)f k (r) 



,i(u>k— WL)(t-t') 



(IV.24) 



where we have omitted the step function since it auto- 
matically gives unity for the integration limits we are 
using here. We will now relate this infinitely short prop- 
agator to some already known parameter. If we go back 
and consider the general result for the equal-space com- 
mutator, this may in terms of the basis- functions {f k } 
be written as: 



[D<->(r,t);D<+>(r,f)] = 
2 



X;«kS(r)&(r) 



3 !(u k -U L )(t-t') 



(IV.25) 
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Comparing with pV.24[) we immediately get a formal 
relationship between this commutator and the infinitely 
short propagator 

(^7-&*.) [D(-)(r,f);D(+)(r,0] 

= -fie i^c 2 P ( ->(r,t|r,t / )- 
(IV.26) 

Using Eq. (|II.32j> this relation can also be written as 



p(-\r,t\r,t') = ±(±-iu> h )fj* t (r,t,t'). 



2c2Ut , dV.27) 

To illustrate how the indefinitely short propagator en- 
ters into the equations we will again consider the sec- 
ond order term in the spin equation represented in Eq. 
(jrV.16j) . The term prior to spatial average is given as 

rt' 



i(3c\c 
he 



dt'J x 



to 



dt 



»d 3 r'{p(-\r,ty,t"Y 

^[J]*D«(r',0} xf(r/) 

(IV.28) 

After spatial average we get two terms, representing the 
coherent and the incoherent interaction. The incoherent 
interaction may then be written as 



^ [ dt' [ dt" J(r) x 

2^o J t0 J t0 



\dt" 



r(r,i',<")V*[J]^ '(r,t")\ xU^M') 

(IV.29) 

To simplify notation, we have signified spatial averaging 
with calligraphic letters, e.g. (D(r, i)) s . a . = 2?(r, i). This 
convention will be used in the remainder of this article. 

We have now developed all the necessary theoretical 
tools to describe the system. In the next section we shall 
use these tools to discuss a pertubative expansion of the 
evolution of the system. 



V. TIME EVOLUTION 

This section is divided into three parts. In the first part 
we examine the general behaviour of the atomic spin in 
the presence of a light field. The aim is to understand the 
effect of the loops introduced in the Feynman diagrams. 
In the second part we consider the light field and we show 
how the theory introduce a decay of the field strength of 
the light as it interacts with the atoms. Again this is con- 
nected to the loops introduced in the Feynman diagrams. 
Finally we will introduce and discuss Stokes operators, 
which are the appropriate operators for describing the 
experiments in Ref . [1, [H, [l(| • 

A. Evolution of the spin 

In this section we will consider the spin equation in de- 
tail for the simple interaction (IIV.10|) . We will begin our 



analysis by considering the first order term in the per- 
turbative expansion of the solution to the spin equation, 
formally given by the diagram 



(V.l) 



This term gives no extra contributions when doing the 
spatial averaging, and we readily write down the expres- 
sion describing this term 



%MV) I' 'dt' J(r,i ) x (-b { -\r,t>) x ©S T '(r,f; 



he 



(V.2) 



We now continue with the second order terms repre- 
sented by the following Feynman diagrams 






(V.3) 



When taking spatial average of these terms, we have 
argued that the first two diagrams will give an addi- 
tional set of Feynman diagrams containing loops and 
stars. It still remains to consider the last diagram of 
Fig (|V.3|) . representing two photons interacting with 
the same atom at time t and t' . In this diagram it 
is necessary to pay special attention to the case where 
the two interactions happen at the same time t — 
t' . The contribution of this term is proportional to 
V { ~ ] {t")T> {+) (t")T> ( - ) (t')V {+) (f) which is not normal- 
ordered, and it will be convenient to separate it into 
normal-ordered terms. When commuting £>(~) (r, t") and 
X> (+) (r,i') we once again get an infinitely short propa- 
gator c.f. (|II.32[) . This extra term we will denote by a 
filled star with a loop. This commutator term will pro- 
duce an interaction which is linear in the field intensity 
(involves X>^~^X>^ + - ) ) whereas the normally ordered term 
(D(-)d(-)dWdW) will be quadratic in the intensity. 
Ignoring for now this quadratic term as well as the coher- 
ent interactions, the second order diagrams for the spin 
equation after spatial average reads 
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(V.4) 



(V.5) 



The interpretation of the diagrams is given below. 

To simplify the expression we will make the slowly 
varying envelope approximation which simplifies Eq. 
jE27| to 



^-)(r,t|r,0 



2c 2 



-v 



f (r,M') 



(V.6) 



Secondly we shall evaluate 77 in a local density approx- 
imation, where we assume that rj(r,t,t') is the same as 
if we were in an infinite medium with a constant density 
p(r) and spin density J(r). By doing this we ignore the 
reflection of the field on the surface of the ensemble or 
other inhomogcncitics. The infinitely short propagator 
which expresses the amplitude for the field to be found 
at the same position at some later time, therefore be- 
comes a delta-function in time. This approximation is 
valid provided that the diffraction matrix M.(v) varies 
slowly on the scale of the wavelength of the light. Fur- 
thermore f/(r, t, t') also contain the Lamb shift which we 
ignore for simplicity. A detailed calculation of fj is pre- 
sented in Appendix [B] where we find 



p(-)(r,i|r,0 



-iS(t - 1') 



q\\(t) 
g±(r) -ig r {r) 
igr(r) g±(r) 



2fc^K-) (r) , 



(V.7) 



where the coefficients Qu, g± and gr may be found in Eq. 
(|B.13I) . Here the result is given in an Euclidean basis, 
where J is assumed to be along the the x-axis. The result 
may also be expressed in a coordinate-independent form 
as 



-iS(t-t') r 4 
- 2 [g±{r) - «7(r)jx 

+ [eii(r)-e±(r)]i(j-}. 

(V.8) 



where j is a unit vector parallel to J. This infinitely short 
propagator is inserted into the second-order terms in the 
spin equation. The second-order incoherent interaction 
given in Eq. (|V.5|) then reads 



He 



dt'l cicdJ^-^-^J • X>[, +) ) - 2><~>(j ■ A^V { + ] ) + H, 



it, 



A^T> { -\3 ■ T> { +) ) - (X><, _) • V { Q +) )A^3 - Tr[A^}V { ~\3 ■ T> { +) ) + T> { ~\j ■ JS^T> { ^ ] ) + H.c. 



(V.9) 



where we have suppressed the space and time dependen- In the simple case, where the matrix is propor- 

cies. tional to the identity matrix, (g-p w 0, Qu ~ g± = g), 
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which is the case to lowest order, the terms proportional 
to C\Cq cancels and the expression reduces to 



2he 



dt' 



>(+) 



(-), 



H.c. 



(V.10) 



This term scale with the power of the incident light, and 
linearly polarized light will affect the spin component 
parallel to the field with twice the rate than the perpen- 
dicular spin components. To see this we may introduce a 
decay-rate Tx>, and writing expression (|V.10[) on a differ- 
ential form, we thus see that the term indeed describes a 
decay of the spin-components. 



dfJ z = — F-r>J z 



(V.lla) 
(V.llb) 
(V.llc) 



and where we have assumed that the light is linearly po- 
larized in the ^-direction. 

Let us now turn to the coherent part of the interac- 
tion represented by the Feynman diagrams in Eq. (|V.3|) • 
The first two terms containing a dot are by construction 
very small, and will vanish when taking the quantum me- 
chanical average, as discussed in Sec. IIVBI The only im- 
portant second-order coherent interaction is therefore the 
following Feynman diagram for normal-ordered fields. 




(V.12) 



where 



_/? 2 c 2 j> (_) {+ ) 
ne 



Suppressing the spatial dependence of the displaced 
electric field, this normal-ordered coherent interaction is 
given in vector representation by 



dt' J T)[ V)(*o \t) ■ V { + \t')) (j • -b { + \t)) (i> l \t) ■ T>[\U)) (Jo • V { + \t))v { + \t>) + H, 



(V.13) 



In the case of linearly polarized light, say T> || e x this 
term vanishes, but this is in general not the case. In Sec. 
IVII we examine the term in some simplified system. 



B. Evolution of the light 

The treatment of the displaced electric field is similar 
to the spin, but there are a few important differences. Let 
us consider the negative-frequency part of the field, and 
write the expansion of the displaced electric field ignoring 
for now the evolution of the spin 




When we take spatial average of diagrams like these, 
we introduce delta- function correlations between vertex 
points. So far we have treated the atoms in the ideal 
gas approximation, where we ignore any correlation in 
the position of the atoms but in reality we should in- 
clude a short-range correlation functions describing that 
two different atoms cannot be at the same position. In 
Appendix [D] we show that including this leads to the 
Lorentz-Lorenz or Clausius-Mossotti relation. In the fol- 
lowing we will only discuss loops, where two consecutive 
vertex points are evaluated for the same atom. Since we 
have subtracted the quantum mechanical average from 
the vertex, no first-order vertex will give a contribu- 
tion to the evolution of the light, and therefore these 
second-order loop diagrams are the most important ef- 
fects apart from the diffraction effects included in the 
mode-functions {f q }. Later in section fV Dl we shall dis- 
cuss the operator nature of the light field and then we 
keep the first-order vertex in the calculations. In the 
current approximation Eq. (|V.14p reduces to 
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3k 




We shall postpone the analysis of this term and discuss it 
in connection with relating the fields to photon counting 
operators below. 



(V.15) 

We have here introduced an interaction denoted by 
a hatched pentagon which scales with (3 2 pkl, and de- 
scribes two connected by the infinitely short propaga- 
tor. Using the results for the infinitely short propagator, 
and taking quantum mechanical average this interaction 
reads on matrix form 



r l|( r ) 

i(3 2 p{r) r ±!l (r) iT r (r) =iM n (r), 




r ±!l (r) iT r (r) 
-*r r (r) r_L )2 (r) 

(V.16) 

where the coefficients entering the matrix are given by 
r,|(r) =c 2 J 4 g n + c 2 g ± (J 2 + J 2 ), (V.17a) 

r_L,i(r) =CoJ 4 g ± + 2c c 1 g r J 2 J x + c 2 (q\\J% + g±J 2 ), 

(V.17b) 

r±, 2 (r) =c%J 4 Q± + 2c oCl g r J 2 J x + dl(Q\\J% + Q±J%), 

(V.17c) 

c 2 

rr(r) =e±2cic 3 2 J x ~ Q\\ y Jx + £>r( c o j2 + cjj 2 ). 

(V.17d) 



We have here suppressed the spatial dependence to 
shorten notation. The series in Eq. (|V.15[) can be in- 
cluded in the differential equation describing the dis- 
placed electric field, 



(2i^j t -Lol + c 2 V x V x [M*(t) + i .M' t (r)])£> M (r,t) 
= c 2 / d 3 r V x V x l T (r,r') • ^[Jj^D^^r',*), 



C. Photon counting and Stokes operators 

So far we have mainly been concerned with calculating 
the field D(r,t). For experiments which eventually in- 
volves counting photons we are more interested in quan- 
tities like photon flux, and in particular the flux in some 
particular polarizational state. We shall now discuss how 
to desribe such photon counting experiments within our 
theory. 

The general idea in this subsection is that we shall 
assume that we are able to measure the light-flux in a 
certain spatial mode by projecting the light field onto the 
mode and then integrating the flux of the light field at 
some detector plane, that we assume to be far away from 
the atomic ensemble. We will formulate such a measuring 
process in terms of an inner product, 

«#r,t)Mr,t)»= / dt f d 2 r^{v,t)-^{v,t). 

J-oc Jm. 2 

(V.20) 

We assume that the fields in general have some axis of 
propagation say i"||. The spatial integral is then per- 
formed in some plane perpendicular to this axis at some 
point rii on this axis. This measuring process could be 
realized by e.g. sending the light field through a single 
mode optical fibre prior to detection. 

We are interested in the polarization of the field which 
is conveniently described by the so called Stokes opera- 
tors defined below. These operators can be derived from 
a Stokes generator defined in a bra-ket-notation by 



g=|D(-)(r,*)))«D(-)(r,*)|, 



(V.21) 



(V.18) which we represent as the following diagram 



where the perturbation is modified accordingly. Because 
of the anti-Hcrmitian matrix, we see that these types 
of loop diagrams correspond to a decay of the field, i.e. 
the differential operator on the left side describes the 
propagation through a lossy medium. On the basis of 
this analysis and the analysis in Sec. IV Al we thus link the 
loops in the Feynman diagrams with the decay associated 
with spontaneous emission. 

It remains to discuss the effect of light interacting with 
an atom that was previously subject to an interaction 
such that the atomic spin state has been changed. In 
terms of Feynman diagrams this is described as 




(V.22) 



Measuring certain light-modes according to the inner 
product in Eq. (|V.20[) . correspond to picking out a cer- 
tain matrix element of the Stokes generator. As an ex- 
ample we assume that in some experiment we are able 
to measure the photon flux of some linear polarization 
in some mode say fq,x 

(r, t) after the interaction with 
the atoms. The time dependence is here fq,a(r, t) = 



The integrated photon flux mea- 



sured at the detector plane, is then given by 



(V.19) 
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Kx\s\i*j 



(V.23) 
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where we normalize the outcome to count the number 
of photons. We have here taken a spatial average of the 
Stokes generator as indicated by the calligraphic font. 

Expanding this operator to second order, gives an ad- 
ditional term not covered by the analysis above. This 
extra term describes a process where both the negative 
frequency part and the positive frequency part of the dis- 
placed electric field interacts with the same atom. This 
extra term comes from the following contribution to the 
Stokes generator 



where K 



2c' 



Using commutation relations for the 




(V.24) 



When taking the spatial average of this term we again 
generate a term representing that the interaction hap- 
pens at the same point. This particular term would not 
have been there if we only considered the spatial average 
of the displaced electric field. The generated term we will 
illustrate as 





s.a. 




(V.25) 



We constructed the interaction represented in the Feyn- 
man diagram as a gray circle, such that when taking the 
quantum mechanical average the term vanish. The new 
term generated when taking the spatial average, given 
as the lower right diagram of Eq. (|V.25|) . describe the 
square of the fluctuations which is not vanishing. This 
was also the case for the terms containing the infinitely 
short propagator. The new term however differs from the 
second order terms containing the infinitely short propa- 
gators because here we need to use the full macroscopic 
propagator. To calculate the effect of this term in detail, 
we therefore need to have an expression for the spatial 
modes describing the system. We will consider this term 
for a simplified system in Sec IV Dl 

To describe the experiments in Ref. [l(| it is conve- 
nient to define a set of polarization dependent photon 
counting operators denoted as Stokes operators. These 
are defined in accordance with Eq. (|V.23j) as 



K 

y 

K 

y 

K 
2? 



;<f*|cS|f*,>: 



;<f*l«s|f*,) 



(V.26a) 
(V.26b) 
(V.26c) 



creation and annihilation operators these Stokes opera- 
tors are seen to have the commutation relations for an- 
gular momentum operators. 



i e. 



q,q 



(V.27) 



We will calculate and discuss these Stokes operators to 
second order in the coupling coefficient (3 in the following. 



D. Calculation of Stokes operators 

In this section we shall calculate the Stokes operators 
to second order. In the experiments in Ref. [1, H, [l(| 
the Stokes operators are measured by sending the light 
onto polarizing beamsplitters followed by a measurement 
of the difference in the intensity of the two outputs. For 
instance if we take the indices q and q' to refer to the x 
and y polarizations of the light, the operator s^' 1 ' in Eq. 
(|V.26j) can be measured by measuring the difference in 
the intensity of the x and y polarizations. The remaining 
operators and Sg' y can respectively be related to the 
difference intensity with the polarizing beam splitter ro- 
tated by 45° and the difference intensity between the two 
circular polarizations. For a general light beam, however, 
diffraction will cause the polarization of the light to de- 
pend on the spatial position and there is no well defined 
polarization. The simple measurement scheme is thus 
only applicable in the paraxial approximation, where we 
can separate out a position independent polarization vec- 
tor. Far away from the ensemble we will therefore assume 
a paraxial approximation. That is, the mode- functions 
f q (r, t) and f q '(r,£) describing the Stokes operators far 
away from the atomic ensemble resemble plane waves 
with transverse profiles that change slowly compared to 
the wavelength. The detector plane is placed far away 
from the atomic ensemble, and at this plane we will as- 
sume that the general set of basis-functions {f q } can be 
approximated as 



f q (r) = J= U n (v ± ) ej e ik * 

V ^7T 



(V.28) 



We have here set the direction of propagation to be 
along the z-axis. The index q are now given as the set 
q = (fc, n, j), where k is some wavenumber, n is an index 
referring to the transverse shape of the mode described 
by the scalar- field U n (r), and j describes the polarization 
of this mode, that can be either x- or y-polarized. The 
completeness relation Eq. (|II.29[) . and orthonormality 
condition in this approximation thus gives 



d 2 r 1 _U* n {r 1 _)U n ,{r 1 _) = 5 m 



(V.29a) 
(V.29b) 



and the dispersion relation Eq. pi.!7a|) at the detector 
plane is w q = c 2 k 2 . 
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Figure 2: Schematic setup. We assume that away from the 
ensemble, the light-mode resembles a plane- wave with some 
transverse profile. A set of lenses focus the beam down into 
the ensemble. 

The paraxial approximation above is convenient for ex- 
pressing the measured observable in terms of the polar- 
ization of the field, but may not be sufficient to accurately 
describe experiments, where tightly focused beams are 
used. We shall therefore only assume this approxima- 
tion to be applicable far away from the sample, and not 
necessarily inside the ensemble. Physically this could cor- 
respond to a situation, where an initially paraxial beam 
is focused onto the ensemble with a lens and converted 
back into a paraxial beam after the interaction by an- 
other lens, as shown in Fig. [21 A similar treatment was 
used in Ref. [3ll |. 

Inside the ensemble we make the much weaker approx- 
imation that the set of spatial mode functions U nq (r) is 
independent of the polarization of the field, so that the 
set f q (r) is given by 

fq(r) = -^=U nq (r)e 3 (r). (V.30) 

The mode U nq (r) now takes into account that the spa- 
tial shape of the beam may change through the ensem- 
ble, and likewise the polarization vector e.,-(r), which we 



shall assume to be real-valued. The index j will still be 
either x or y, corresponding to the polarization of the 
mode far away from the sample, but the vector e.,(r) 
will not necessarily be parallel to the x or the y axis. 
A more general description of the mode-functions would 
include a dependence of the polarization vector ej on the 
polarization state U mq (r), i.e., e m j(r). The correction 
this generalization gives to the Stokes operators, is pre- 
sented in Appendix [H] in relation to Sec. IVI CI When 
we make the relevant calculations to describe the Stokes 
operators defined in Eq. (|V.26[) . we will chose to con- 
sider modes corresponding to the index q = (k, m, x) and 
q' = (k, m', y). We note that the set {f q } defined in this 
way is in general not complete, since, e.g, the assump- 
tion that the polarization vector is independent of the 
transverse mode number applies in the paraxial approxi- 
mation but does not apply in general. When calculating 
the effect on the forward scattered field to first order we 
only get contributions from the near paraxial modes in 
the forward direction. When we go to second order there 
will, however, be effects of all the transverse modes, and 
in this case a correct treatment requires a more accu- 
rate treatment of the complete set of modes. Above we 
have already employed such a more general set of modes, 
when we discussed the effect of spontaneous emission, 
which involve all the transverse modes. In addition to 
this, a more accurate set of modes is also required for 
describing the effect of dipole-dipole interactions, which 
also involves all the transverse mode. 

We will in the following calculate the Stokes operators 
in the limit described above. Diagrams containing a loop, 
we will not discuss, since these only leads to a decay of 
the light which we have discussed earlier. After taking 
spatial average the diagrams in consideration are 




^v^^ 




Let us begin our discussion of this perturbation series 
by considering the first term on the right hand side of 
equation (|V.31[) . This term is the zeroth-order term of 
the Stokes generator S^°\ In the far-field limit z 



(V.31) 



I 

the matrix-element we need to calculate is 



dtd 2 r ± — L U m (r ± ) ej e lkz -^* • 

i V27T 



(V.32) 
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and thus gives us 

K{{f* kmj (r, f)| \r kmrj , (r, t))) = 4 m3 «« j( . (V.33) 

The zeroth order Stokes operator s qq for q = (k,m,x) 
and q' = (A;, m', y) gives 



•;qq' ~mm' /at 



mx O,kmx 



a 



km'y akm 'v 



(V.34a) 



The two remaining zeroth order Stokes operators are 
found accordingly, 



'ylliril 

s 2 



1 



), (V.34b) 

J " m ' = 2i ( a *mx a *m'j, " oL'y a *™») ' ( V - 34c ) 

In the following we will calculate the first-order compo- 
nents of the Stokes operators. We assume the quantum 
mechanical average of the atomic spin J to be parallel 
the x-axis. The relevant interaction matrix can in this 
case be written 



i [J] = ici/3 



J z (r) -J y (r) 
-J*(r) 



Jy(r) 











(V.35) 



and after spatial averaging we simply write 



-in >!'ir< | J y (r) | - 



(V.36) 



The second and the third term on the right hand side of 
Eq. (|V.31[) are the first order terms of the Stokes gen- 
erator, . To calculate the contribution to the Stokes 
operators from these terms we have to evaluate the ex- 
pression 

((r krnj (r, t)\ c 2 j£ dtd 3 r' P(")(r, t\r' , t')- 

^[J]X>S _) (r',0». 

(V.37) 

The initial time to we will set to — oo, and because we 
assume our detector plane to be infinitely far away from 
the atomic ensemble, we can take t — ► oo. Using the 
expression for the set {f q } given by Eq. (|V.28|1 for the 
detector plane and Eq. (|V.30|) inside the ensemble, Eq. 
(|V.37[) reduces to 



nl 

(V.38) 
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where 



e^r') = U km {vrU kn {v')e 3 {v) ■ [ J y (r>) x e,(r')] 

Mr') J 

= *r(r') [SiJ jy - S JX S ly ] [ | J v {r>) | • e 2 (r')] , 

Ur')J 

(V.39) 



with 



^r(r')=U km (r'yU kn (r'). 



(V.40) 



In the final equality we have introduced the local basis 
vector e z (r) = e x (r) x e y (r). The effect of the first-order 
term of the Stokes generator to the Stokes operators 
thus reads 



k LCl (3 J d 3 r'^p(r') i{ejr(r')*«L«^' 

nl 

+ Qf'i n (r')al mj a knl }. (V.41) 



The remaining terms of the right hand side of Eq. 
(|V.3ip . that is the second-order terms, can be calculated 
in a similar way. The results may be found in Appendix 
|U The calculations given in Eq. (|V\4Tj) . (|KT|) , (pO|) and 
(|E.5p is the starting-point for a discussion of the dynam- 
ics of the system subject to a general light field of many 
modes. 

The description that we have used here, where we de- 
fine the Stokes operators in term of expectation value 
between different orthogonal modes, is very convenient 
for a theoretical description of the process. It does, how- 
ever, not directly correspond to the experimentally mea- 
sured observables unless one, e.g., separates out partic- 
ular modes with single mode optical fibers. We shall 
therefore defer the discussion of the consequences of these 
results to the next section, where we use these result to 
calculate the evolution of observables more relevant to 
experiments. 

We will now give the equation for the atomic spin. The 
incoherent terms describing decay due to spontaneous 
emission have already been discussed. Here we will con- 
sider the coherent interaction up to second order in the 
perturbation series. Below we show the diagrammatic 
representation of the coherent perturbation series for the 
atomic spin up to second order. 




We will denote the first order term in the expansion, 
Eq. (fV\42)l . as J (1) . Employing again the approxima- 
tions done in the previous calculations, that is, using the 
set of light modes {f q } given in Eq. (|V.30[) and setting 
the initial time to — oo and the final time to oo, the term 
can be written 

= -/? Cl fc L J2 ^ m 'W(j(r) x e ,(r)) 

kmm' 

U.j , _„ t , 1 

2j [ a krnx akm 'y a kmy akm ' x 

= -fJ Cl k L (j(r) x e z (r)) 

kmm' 

|Rc[*^ m '(r)]s™ m ' + Im^'^V)]^™ 1 "'}- 

(V.43) 

We notice that compared to the simple theory in Ref. 
[2^ | there is an additional term proportional to the imag- 
inary part of the function ij/ m ™ 1 (r) . A similar correction 
can also be found for the Stokes operators for the light. 
Also notice that the dynamics of the spin to first order 
happens in a plane orthogonal to the vector e z (r). This 
is the reason why the term in Eq. (|E.3jl vanish, since 
there we are considering the effect of the dynamics of the 
atomic spin on an axis parallel to the e z (r)-vector. The 
calculation of the second-order terms is presented in Ap- 
pendix |FJ In the following section we will examine the 
effect of these calculations under conditions attainable in 
experiments. 

VI. EXPERIMENTAL APPLICATION AND 
VALIDITY 

In this section we shall consider different limits where 
we can reduce our general theory to a theory resembling 
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(V.42) 



I 

the simple description obtained in one dimensional the- 
ories 0, H^l • Furthermore we discuss the validity of the 
approximations made to arrive at these simple limits as 
well as the validity of our perturbative treatment of the 
interaction. 



A. Measurement procedure 



In the previous section we discussed how our theory 
could be used to calculate Stokes operators correspond- 
ing to specific transverse modes of the field. While such a 
treatment is appealing from a theoretically perspective, 
it is less desirable experimentally, since the isolation of 
single transverse modes is complicated (although it could 
be done by passing the light through single mode opti- 
cal fibers). Here we shall therefore express our result 
in terms of a simpler experimental procedure. Suppose 
that the detections is performed by sending the light onto 
a polarizing beamsplitter and recording the intensity of 
the two output port with two cameras. The difference 
between the intensities can now be used to define posi- 
tion dependent Stokes operators §i(r±), i.e., si(rj_) cor- 
responds to the difference in intensity between x and y 
polarization at position rj_ in the detector plane. Sim- 
ilarly S2(r_i_) and S3( r -iJ can, respectively, be related to 
the difference intensity with the polarizer rotated by 45° 
and the difference intensity between the two circular po- 
larizations. These operators may in general be deter- 
mined by 
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kmm' 
kmm' 

h(r±) = ^(Um( r ±)at mx akm>yU m >(r±) - U m {r ± )ai my a km , x U m/ {r ± )^ 

kmm' 



Below we shall derive expressions for the operators 
(jVI.ip and discuss how to implement a light-matter quan- 
tum interface based on these operators. In subsec. IVI Bl 
we for simplicity first consider an extreme paraxial limit, 
where we assume that essentially no diffraction occurs 
during the propagation. In this limit the dynamics be- 
comes extremely simple. In subsec. IVI CI we consider 
a more interesting limit, where we may have multiple 
modes which may experience diffraction. Here we show 
that measurement of the operators Si(r±) still allows us 
to simplify the dynamics of the system. In a suitable 
limit we find a simple two mode transformation between 
transverse modes of the light field and single modes of 
the atomic ensembles. 



B. Extreme paraxial approximation 

In the extreme paraxial approximation, we completely 
ignore any dynamics transverse to the propagation direc- 

I 



h,out(rx) 



S2,out{r±) 



h,out{r±) 



In this limit we see that the Stokes operator S3 is de- 
coupled from the coherent dynamics of the system, and 
only evolves due to spontaneous emission [derived in Eq. 



(VI. la) 
(VI. lb) 
(VI.lc) 



I 

tion of the light modes and approximate the set of modes 
{f q } with Eq. (|V.28[) throughout the ensemble. Since the 
typical distance for diffraction is given by Id ~ A/X, the 
condition for the validity of this approximation is L <C Id, 
or expressed in terms of the Fresnel number T 1. 



The full expressions for the Stokes operators are quite 
involved, and we therefore leave the incoherent part of 
the evolution to Appendix [G] Keeping only the coherent 
part of the interaction, we find the Stokes operators to 
second order in the interaction to be 



(VI. 2a) 

(VL2b) 
(VI.2c) 

I 

Similarly we may find the coherent dynamics of the 
atomic spin. Leaving again the incoherent part to Ap- 
pendix [G] we find 



=si,in(rx) - k h ci(3 I dz' p{z' ,vx)J z (z' ,T±)h,in{*±) 
1 



(fc L /? Cl ) 2 WVf/.rxW/.rxl^/.rxJJ.f/.rj^^rx), 



=S2,in(rx) + k h ci(3 I dz'p(z' ,r±)J z (z' ,rj.)Ji,j„(rj_) 
1 



(fc L /? Cl ) 2 // dz'dz"p(z',T ± )p{z",T ± )J z {z',r ± )J z (z",r ± )s 2 ,in(r±), 



= S3,m(r_l_). 
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Jx,out( r ) — Jx,in( r ) ~ /3ci/c L S Jj/,in( r )^3,m( r J 

fc 

■4,o«t(r) =J 3/ ,i„(r) + /3cifc L J a) fa(r)s| iw (rj 



i(/?Cifc L ) 2 51 ^,m( r )s3,m(rx)C»( r ^) 

i(/3cifc L ) 2 ^ ^,m(r)s| )in (r ± )s|' in (r ± ) 



A- A' 



(VI.3a) 

(VI.3b) 
(VI.3c) 
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Analogous to what we found for S3, we see that the op- 
erator J z is decoupled from the coherent dynamics of the 
system. This result can directly be associated to the con- 
servation of angular momentum along the z-axis. In the 
extreme paraxial approximation this is true to all orders 
in the coherent interaction. 

The results in Eq. (fVL2|) and (|VL3l) is essentially 
equivalent to the simplified one-dimensional description 
of the system given in Refs. @, . The only difference 
is that the expressions derived here now apply for each 
value of rx, whereas the previous treatments assumed 
the system was transversely homogeneous and only con- 
sidered the variables integrated over r^- 

A further simplification of Eq. (|VI.3[) can be obtained 
if we introduce the rotation vector 



n = /3ci/c L ^^ m (r ± )e^ 



(VI.4) 



With this definition we find that Eq. (|VI.3[) describes 
nothing but a rotation of the spin around the e z -axis 



3 out = Jm + Jin X ft + - (Jj„ X fl) X ft. (VI. 5) 



C. Multi-mode coupling 

In the previous subsection we basically ignored all the 
dynamics transverse to the propagation direction. Now 
we turn to a more interesting situation, where we may 
describe effects associated with diffraction of the light 
beams. Our goal in this section is to find a set of condi- 
tions under which we can have a simple dynamics, where 
the individual transverse modes of the light field talks to 
a single mode of the atomic ensemble. Such an interac- 
tion would enable the storage of information from several 
light modes into spatial modes of the ensemble, e.g., us- 
ing the protocol in Q . The realization of this interaction 
would thus expand the information storage capacity of 
the atomic ensembles. A similar problem is considered 
in Ref. In related work such storage of multimode 

memory has recently been achieved in atomic ensembles 
using electromagnetically induced transparency [33| . 

To achieve simple results in the end, we will here con- 
sider a situation, where we have a strong classical beam 
polarized in the x-direction in a single transverse mode 
U k(r) (denoted by the index o). For the y-polarization 
we, however, include a complete set of modes, which may 



or may not include a term identical to the mode of the 
x-polarization. For the strong mode we will approximate 
a\ ox — dkox — y/ N° ^> 1 where N° is the number of pho- 
tons in this particular mode. Since the Stokes operators 
are dominated by the terms involving the classical com- 
ponent, the only important contributions in the Stokes 
operator (|VI.1|) are the terms containing the strong clas- 
sical mode. Eq. (|VI.1[) are thus approximated by 



4 m) (r±)«2|CWr ± )| 2 7V°, 



(VI.6a) 
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Re[U* ok {rx)U mk {v x ))X^ 

-Im[C/ * fe (r ± )C/ m)fe (r x )]P^ 

(VL6b) 
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+ Im[U: k (v ± )U mk (v ± )]X?), 
(VI.6c) 



where 



Ip =-= 



V2 



y^kmy + ®kmy ) 



P2 1 =■ 
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(VI.7a) 
(VI.7b) 



In order to obtain simple result in the measurement 
process, let us assume that we can choose the mode func- 
tions t/ m fc(r) to be real in the detection plane. This could, 
e.g., be achieved by sending the light through a lens 
which converts the incoming modes into extreme paraxial 
beams as shown in Fig. [5] (note that since we only make 
this assumption in the detection plane, this assumption 
does not restrict the shape inside the ensemble). Experi- 
mentally the operators X m and P ra defined here can then 
be measured by simply integrating the measured §i(r±) 
with a suitable weight function, e.g., 




s 2 (v±)=Xl 



(VI.8) 



where we have used the expansion in (jVI.lj) as well as the 
orthogonality relation of the transverse mode functions 
(IV.29bl) . 
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In our equations of motions we for simplicity only keep terms. The equations of motion for the Stokes operators 
terms to first order in (3 and \/N°, and neglect all other give in this limit 

I 



=-v;;; +/.,. j dV^o ( J y (r) \ ■ c.jdr.- #ru-) 



P%t=PZ + k,Pc A r-f /dVp(r') | J» | ■ e,(r)Im[*r(r)], 
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(VI.9a) 
(VI. 9b) 



where ^ mo is defined in terms of the mode functions U m 
in Eq. ()V.40|) . Employing the same set of approximations 



J 



for the spin equation we find 



Jo«t(r) « J ln (r) + fc L /3ci 




Re[*r(r)]^ - M*3T (j«(r) x e z (r)) . (VI.10) 



r 



Note, that the expressions we have derived here, allow 
for a general set of transverse modes which may experi- 
ence diffraction, and thus go beyond the extreme parax- 
ial approximation made in the previous section. In the 
expressions above we do, however, still use the paraxial 
approximation in Eq. (|V.30[) , where we ignore the depen- 
dence of the polarization vector on the mode number. In 
Appendix [H] we relax this approximation. 

The expressions in Eqs. (|VI.9[) and (jVI.lOp differ from 
the simple results of the last section because of the ex- 
tra terms proportional to Im[ , 5™ (r)]. These terms com- 
plicate the dynamics and, e.g., means that one cannot 
use the protocol in Ref. [9j to store information in the 
ensemble. There are, however, certain limits where the 
extra terms in Eq. (|VI. 10[) disappear. One situation is 
when the mode we are considering in the y-polarization is 
identical to the classical mode in the re-polarization (ex- 
cept from the different orientation of the polarization). 
This situation corresponds to the experimental situation, 
where the weight factor U m /U in Eq. (|VI.8|1 is unity, 
such that the final result is obtained by integrating the 
intensity over the transverse plane. This case therefore 
corresponds the experimental situation where the light is 



detected by photo detectors instead of cameras. In this 
case Im[ , I'™ (r)] vanish identically and the evolution of 
the light operators again resemble the result of the last 
section, where, e.g., the $2 component was conserved, 
which translates into P™ t = P™. Note, however, that 
unlike the situation considered below, the atomic opera- 
tors in this situation gets an admixture of several different 
input light modes, and will not in general reduce to the 
dynamics considered in Ref. @. 

Let us now consider a different limit ideally suited for 
a multi-mode memory. We assume that we are in the 
paraxial approximation, where we can ignore the spatial 
dependence of the polarization vectors. For simplicity we 
also assume that the classical mode U (r) has a uniform 
intensity and that the density is constant over the re- 
gion, where U m is non-zero in the atomic ensemble. We 
furthermore assume that the macroscopic polarization is 
constant and along the x-axis, J x , and finally we assume 
that iff 1710 is real (for a discussion of the validity of this ap- 
proximation we refer to the next subsection) . In the spin 
equation (jVI.lOp we will only keep terms proportional to 
the macroscopic spin component J x . In this situation the 
relevant equations reads 
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(Villa) 
(VI. lib) 

(VI.llc) 

(VI. lid) 



Here the factor exp(— ikz) comes from the classical field 
and cancels the exp(ifcz) dependence of the mode func- 
tion U mi since C/ m exp(— ikz) should be real according to 
the assumption of '5 being real. This set of equations 
can be symmetrized and simplified by introducing a set 
of collective operators 



X' 



—ikz 



P A=\Iyl I d3rJ *(r)U m (r)e- lk ~ 



(VL12a) 
(VI. 12b) 



where L is the length of the ensemble. The coefficients 
here are chosen such that the operators X™ and P™ 
fulfil the standard commutation relation for position and 
momentum 



[X^-iPa ] — i&mm' ■ 

With these definitions Eqs. (jVI.lip reduce to 



^P.out ~^P,in 



HP 



A.in i 



P.out 



P.in i 



x 



A. out 



A, out 



--X 



A,in 



A.in i 



kP, 



P.in ■ 



where 



K =/fi L /?Ci[/o 



N°pJ x L 



(VI. 13) 

(VI.14a) 
(VI. 14b) 
(VI. 14c) 
(VI.14d) 

(VI.14e) 



These equations describe a system where one transverse 
light- mode couples to a single mode of the atomic ensem- 
ble, which in term couple back to the same light mode. 
This two-mode mode dynamics is exactly identical to the 
dynamics derived in Ref. 7| for a single transverse mode. 
The dynamics can thus, e.g., be used to realize a multi- 
mode version of the memory protocol implemented Ref. 
Q. In this protocol P™ ira is stored in the atomic mode 

-V™ out , while at the same time the atomic mode P™ in is 

transferred to the light-mode Xp out , as described by Eq. 

(|VI14[) . After detection of the light operator Xp out one 
can then realize a quantum memory by feeding back the 
measurement result to the atoms as it was shown in Ref. 



D. Validity 

1. Validity of the simple multi-mode dynamics 

In the previous subsection we derived a simple multi- 
mode dynamics useful for making a multi-mode light 
matter quantum interface. For experimental implemen- 
tation of these idea an important question is the validity 
of the approximations leading to Eq. (|VI. 14|) . First of all 
we need that the imaginary part of vI/ om (r) in Eq. (|VI.9|) 
should vanish. Furthermore, in order to define orthog- 
onal spin-modes that do not couple different transverse 
modes, we need |J7 (r)| to be uniform. Taking the classi- 
cal mode to be given by U {v) = U e tkz , where U is real, 
we also need the quantum mode U m (r) to be real- valued 
apart from the e %kz dependence. Let us now take the 
modes U m (r) to be Hermite-Gaussian beams 34]. Such 
modes can be represented by 

tWr) =^H m ( V2^-) H n ( V2^) 

w(z) \ w(z) J \ w(z) J 

^, i[kz— (m+n+1) tanhz/zo] 

x el k(x 2 +y 2 )/2B.(z) e -(x 2 +y 2 )/w 2 (z) 



where 



w(z) =w Q ^Jl + z 2 /z%, 



R(z) =z + ^, 

z 

_nw 2 
*0-— 



(VI.15a) 

(VI. 15b) 
(VI. 15c) 
(VI.15d) 



Here wq is the minimum waist of the beam, k is the wave- 
number, A is the wavelength, B G R is a normalization 
coefficient, and H n is the set of Hermite polynomials. 
The condition that U mn (r) must be real- valued gives the 
conditions 

XR(z) > w 2 (z) \{l + m + n) — | < 1 (VI.16) 

These are in fact equivalent conditions, and introducing 
the Fresnel number T = w 2 {z)/\L we find the condition 



T ^> 1 + m + n. 



(VI.17) 
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2. Validity of perturbation theory 

The theory we have developed in this paper is based on 
perturbation theory in the interaction between light and 
atoms. In this subsection we discuss the limits of valid- 
ity of this perturbative treatment. We will be considering 
worst case scenarios to find the limit, where our perturba- 
tion scries Eq. (|V3T|) and (fQ2|) converge. An impor- 
tant parameter for these estimates will be the effective 
coupling constant for the collective operators n defined 
in Eq. (|VI.14[) . For applications to light-matter quan- 
tum interfaces this parameter should be of order unity. 
As we shall see below, this is still possible without vio- 
lating the applicability of perturbation theory. Another 
important parameter is the optical depth, OD, defined 
by OD ~ p\ 2 L. The optical depth plays an important 
factor when describing the effect of the incoherent inter- 
action, e.g., the spontaneous emission. 

Throughout this work, we have assumed that the 
atomic ensemble is polarized along the rr-axis, so that the 
atomic spin components pJ y , pJ z only carries quantum 
noise. Also we have assumed that the classical compo- 
nent of the light is linearly polarized so that, e.g., circu- 
lar components are governed by quantum noise. These 
assumptions will be important for estimating the terms 
below. 

We first consider the expansion of the light field (|V.31[) , 
and in particular the coherent part of the interaction. 
The effective perturbation coefficient for the first order 
term is found to scale at most as (/3fc L V 'N a) / 'A ~ n/^/Np 
(may be found by estimating Eq. (|V.4ip ). Here A is the 
transverse area of the atomic ensemble, and Np is the 
total number of photons in a pulse. Going to second or- 
der an important term is described in Eq. (|E.3|) . Since 
we are not including the the time evolution of the macro- 
scopic polarization in the average interaction, this term 
has a potential scaling as large as k 2 . We showed, how- 
ever, that in the paraxial approximation the term van- 
ish. Going beyond the paraxial approximation as done in 
Appendix [H] we find that for linearly polarized light the 
scaling is n 2 /y/Np. The last contribution to Eq. (1V.31|) 
is the incoherent interaction considered in Appendix [G] 
The scaling of this effect k 2 ■ (N A /N P )/OD . 

Now we consider the spin series (|V.42j) for a single 
atom. The incoherent part of the evolution of the spin is 
described in Eq. (|V.9[) . and scales as k 2 /OD, it can be ig- 
nored for sufficiently large OD. The first order term scale 
as k/^Na for linearly polarized light. To increase this 
coefficient we need circularly polarized light, which makes 
it interesting to examine the second order term describing 
the change of the polarization of the due to the interac- 
tions with atoms. This process is described in Eq. (jF.ljl . 
which represent the optically induced dipole-dipole inter- 
action. This particular term vanish when we take quan- 
tum mechanical averages, because we have subtracted the 
only non- vanishing component, but we can still calculate 
the root mean square contribution. The effect can then 
be separated into a short range part and a long-range 



contribution. The long range contribution can be esti- 
mated to give a contribution of order n 2 y/d/ (L ■ OD), 
where d is the smallest dimension of the setup, i.e., the 
smaller of the length and the transverse sizes of the beam 
and the ensemble. The short range part actually diverges 
within our present approximations. If, however, we reg- 
ularize the integral by excluding the volume, where the 
dipole-dipole interaction of an excited and a ground state 
atom V ~ 7 A 3 /r 3 is of the same order as the detuning A, 
we find a contribution k 2 'y/ 'A/7^/ 'A/ '(L ■ OD). The justi- 
fication for this regularization is that when we made the 
adiabatic elimination we assumed a constant detuning 
A. This approximation breaks down when two atoms 
are sufficiently close that the dipole-dipole interaction 
is the strongest effect in the problem, in which case it 
is more appropriate to describe the atoms in terms of 
molecular states. Both the short and long range part of 
the interaction are thus small for sufficiently large opti- 
cal depth OD and for sufficiently long ensembles (large 
L). It should, however, be noted that here we have only 
performed a very rough treatment of the dipole-dipole 
interaction, and it would be desirable to make a more 
accurate treatment of the effects of these terms. Also it 
should be noted that the estimates we have performed 
here apply to non-moving atoms, i.e., cold atoms. If we 
include the motion of the atoms, i.e., warm atoms as in 
Refs. 

there will be a reduction of these terms 
because the sign of the interaction will change in time. 

In summary, sufficient requirements for the conver- 
gence of the series for the light fields are 



«1, 



«1, 



K 

OD 



Na 
N P 



< 1, (VI.18) 



and for the spin equation sufficient requirements are 




«1, 



(VI. 19) 



By having many atoms and photons as well as a large 
optical depth, it is thus possible to achieve k *~ 1 without 
violating the applicability of perturbation theory. 

The main idea in this work is to develop a perturbation 
series, where we explicitly take into account the reshaping 
of the light modes caused by the mean effect of the inter- 
action. Let us for comparison compare with the series, 
if the mean effect of the interaction had not been sub- 
tracted. For the Stokes operators the perturbative series 
is given in Eq. (|V.31|) . If we do not subtract the average 
effect of the interaction, the scalar part of the interac- 
tion [the Co component in Eq. (|II.4[) ] will give first order 
corrections to the field of order KyJ Na/Np times the in- 
coming field. With Na ~ Np as it is suggested in Ref. 
0, this term will give a factor of order unity for k ~ 1, 
and this therefore cannot be considered a small term. For 
the calculation of the Stokes operators, however, the two 
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large components in the first order terms in Eq. (|V.31|) 
cancel out. The calculation may thus yield reasonable 
result even without performing the more involved pro- 
cedures described in this article, but the validity of the 
procedure would be questionable. (Some experiments ac- 
tually uses N p ^> Na [H, where this problem may be of 
minor concern). Furthermore, one of the major limiting 
factors identified above, is the dipole-dipole interactions. 
The effect of this term is much more complicated to eval- 
uate if we had not subtracted the average interaction, 
but the term certainly will be larger, because the inter- 
actions in Eq. (|F.1[) would include a non- vanishing term, 
and not just the quantum fluctuations. Again this term 
would thus seriously question the applicability of pertur- 
bation theory. In contrast the present approach allows 
us to rigorously apply perturbation theory in experimen- 
tally relevant regimes. 



VII. CONCLUSION 

In quantum optics the propagation of light through an 
atomic medium is often described in a one-dimensional 
approximation, where one completely ignores the trans- 
verse structure of the beam and only considers the longi- 
tudinal propagation. In this paper we have investigated 
the validity of this approximation by developing a full 
three-dimensional theory describing the interaction. The 
challenge in this work has been to develop a theory ca- 
pable of describing the microscopic interaction with a 
single atoms as well as macroscopic effects such as the 
diffraction of the laser beam caused by the refractive in- 
dex of the gas. In essence the theory we have developed 
here includes both the micro- and macroscopic effect by 
separating the interaction into an average part and the 
fluctuation from the average. In this formulation macro- 
scopic effects such as diffraction are naturally associated 
with the average part whereas the microscopic fluctua- 
tions describe processes such as the mapping of quantum 
fluctuations between light and atoms. Furthermore we 
have shown that spontaneous emission from the atoms 
naturally appear as an effect caused by the fluctuations 
associated with the point particle nature and the random 
positions of the atoms. 

Based on our separation into the average and the fluc- 
tuations we have developed a perturbative expansion in 
the fluctuations. The advantage of this procedure is that 
it has a wider region of applicability than a direct pertur- 
bative treatment. For instance in an experimental setup 
an index of refraction of the gas just change of the beam 
profile which often only has a minor effect on the experi- 
ment. On the other hand, such 'trivial' effects may have 
a large influence on the theoretical calculation. If one 
considers perturbation theory based on the vacuum so- 
lutions to the wave equation, the perturbative expansion 
will include all the terms responsible for the reshaping 
of the beam, and this may break the validity of pertur- 
bation theory. On the other hand our theory performs 



perturbation theory on modes which are solutions to the 
wave equation including the index of refraction of the gas. 
Our theory is thus applicable even for situations where 
the beam is considerably distorted by the refractive index 
of the gas. 

A major motivation for this work has been to investi- 
gate the validity of the one-dimensional approximation in 
the description of the experiments in Refs. [81 \9L lid. 
In Sec. IVII we explicitly considered some situations where 
we could reduce our general theory to a theory resembling 
the one used to describe these exp eriments in the one di- 
mensional approximation 0, l29l |. To achieve a simple 
description resembling the previous theories, an essen- 
tial requirement is that we are in the paraxial approxi- 
mation. If we are not in this limit, the polarization of 
the light change as its propagate through the ensemble, 
which complicates the interaction with the atoms. Fur- 
thermore, for the particular interaction considered here, 
we also find it to be desirable to be in a regime where 
the Fresnel number is much larger than unity T 3> 1. In 
these limits our theory essentially reproduce the results of 
the simple theory. The only difference is that instead of 
the vacuum mode functions, the mode functions appear- 
ing in the theory should represent the modes, which are 
solutions to the diffraction problem including the index 
of refraction of the gas. 

In the present paper we have mainly focused on de- 
veloping the theory and deriving how the usual approxi- 
mations arise from our more complicated approach. The 
theory is, however, fully consistent and thus capable of 
including any higher order corrections not previously in- 
cluded in the theoretical description. In particular it 
could be interesting to study the effect of light induced 
dipole-dipole interactions. While such processes may not 
be relevant for understanding the current experiments, 
they may play an important role in future experiments, 
e.g., with Bose-Einstein condensates, where the density 
may be fairly high. Another interesting extension of our 
theory could be to study different types of interactions 
such as for instance electromagnetically induced trans- 
parency [n|. 
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Appendix A: ADIABATIC ELIMINATION 

In this appendix we derive an effective Hamiltonian in- 
volving only the atomic ground state. The Hamiltonian 
can be expanded on the complete set of states de- 
scribing the atom. Let such a set be comprised of a set 
of exited states {|ej)} and a set of ground states {|<?i)} 
so that the Hamiltonian reads 

ft = X)( w J' + w °)l e ^^l +12 u o\9i)(9i\ +Wmt. (A.l) 



For convenience we have here set h = 1 and only consider 
a single atom. The set of ground states are assumed to 
have the same energy, ujq and uij is the transition fre- 
quency from the ground state to the exited state | e j- ) . 
The interaction Hamiltonian is given in Eq. (fTTT3|) . and 
when expanded on the set of internal atomic states it 
reads 



n^ = --Yp { -\t)-{9i\P\e j )\g i ){e 



+ {e j \P\g i )\e j ){g i \-i>^(t), (A.2) 

where we have used the rotating wave approximation as 
well as the fact that the matrix elements (e,-|P|ej<) and 
{gi\P\gi') vanish. To shorten the notation we suppress 
the spatial dependence. We will use that the displaced 
electric field primarily oscillate at the laser frequency, 
and change to the interaction picture 



6(->(t) a < 



(A.3) 



Using Heisenberg's equations of motion we may derive 
an equation of motion for |<7j)(ej| 



^lffi)( e jl = ~ iA j\9t)( e j\ ~ ^-^2{(ej\P\9i)\ef)(e 

j' 

-(e J \P\g l ,)\g l )(g l/ \}-& + \t), 



(A.4) 

where is slowly varying. In the limit of weak driving 
we may set J^|<7i)(ej| = 0, and obtain an approximate 
solution 

\9i){ej\ * ^J^{e j \P\g i ,)\g i ){g i ,\-i)^(t), (A.5) 

i' 

where we have neglected the exited state population. The 
atomic part of the Hamiltonian can be written 

fto =^2A j \e j ){g \go){e j \ +J2uj L \g l ){g t \ 

J i 

+ E( W ° ~ ^)(|e i )(e J | + \gi){gi\), 

(A.6) 



where | <7o) is any ground state. By inserting expression 
(|A.5jl and the Hermitian conjugate into Eq. (|A.2|) and 
(|A.6jl we find the simple result 



e V ^7 £ 0^7 



H 



| 9i )(^|(( ej |P|^)-DW(t)). (A.7) 

(neglecting a zero-point energy term in the Hamiltonian) . 
We may now identify the matrix operator V[J] 

= E ^-(^|P| e j)( e i|P|^)l^)(^l (A.8) 

jii' ^ 

and we immediately get the result stated in equation 
(|II.3|> . The notation "•" in this expression means usual 
vector product with the vector to the right. Furthermore 
we may also find the relation between the polarization 
and the displaced electric field 



p(-)(f)=^|e i )(e i |P| 5 ,)(^| 

ij 

jii' "* 

= F*[J]D<->( t ). (A.9) 

We have here only written the positively oscillating com- 
ponent, the negatively oscillating component is found by 
Hermitian conjugation, which from equation (IA.8j) is the 
same as transposition of the matrix. 



Appendix B: CALCULATION OF INFINITELY 
SHORT PROPAGATOR 

In this appendix we calculate the infinitely short prop- 
agator in the local density approximation. We will for 
simplicity only consider the simple interaction given by 

V[J] = /3p(r)( Co J(r) 2 - j Cl J(r) x ) . (B.l) 

We further shorten the notation by introducing the coef- 
ficients a = 1 — /3p(r)c J(r) 2 and a\ = /3p(r)ci|J(r)|. 

If we Fourier-transform equation (|II.17|) . the equation 
we wish to solve is 



k x k x (n I iaijx)e k = — 



2 k 2 



(B.2a) 
(B.2b) 



where the vectors k and j are unit vectors representing re- 
spectively the direction of the plane wave solution and the 
orientation of the atomic spin. The solutions to the above 
equations is the following set of polarization-vectors 



/ i x k kx(ixk)\ . , N 
V|jxk| |kx(jxk)K ±V ' 



(B.3) 
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where vi and V2 are unit vectors given by the first and 
second fraction respectively. The normalization constant 
N± is determined by using the inner product in Eq. 
pi.!8[) . In this way we find the real space representa- 
tion of the basis- functions fk(r) 



f£(r) 



1 



(vi±«r 2 )e <k,p . (B.4) 



/ 2(2^)3(a ±ai(j-k)) 
The dispersion relation is then derived from (|B.2a[) 

c 2 fc 2 (a ±ai(j-k)). (B.5) 



w 2 ± 



The infinitely short propagator can then be calculated 
to be the following 



2lu,c 2 



se{+,-} 



]T / d 3 fc^ s (f s k (r))*f s k (r). 



(B.6) 



We introduce the matrix given by the following juxtapo- 
sition: 



M(k,j,s) = (vi - isv 2 )(vi +isv 2 ) 



(B.7) 



Changing to spherical coordinates and making the sub- 
stitutions x — cos 9 and k' = ky/1 — ao + sa±x as well as 
using the dispersion relations given in equation (|B.5[) the 
integral reduce to 



P ( -\r,t-t') 



2lu, c 2 



E 



«e{+,-} 



OO />1 />27T 

dk' dx 1 

n 7-1 JO 



2U4 



C 2 k 



2(27r) 3 (a + sa lX fl 2 
I 



Mix^^y^'^' 2 -^/^. (B. 



Neglecting the denpendence of k! outside the exponential 
and using that the difference k' 2 — k 2 for large A: L runs 
from —00 to 00, the k' integral gives a delta- function 
in time. Including the <f> integration in a matrix M we 
finally get 



-ikl5{t-t') 
16^V 



»e{+,-} 



1 dx m ^ s K 

1 (ao + saix) 2 



with the matrix M given by 

M(X, S) = 7T 



2(1 -x 2 ) 
1 + x 2 2isx 
-2isx 1 + x 2 



(B.10) 



The s-sum is evaluated by substitution in the integral and 
the final expression for the infinitely short propagator is 

= , -ik?5(t-tf) f 1 , M(x 
P(->(r,t-t') = kA-z -/ dx 



8ttc 2 



_l Tr(a + aix) 5 / 2 

(B.ll) 



These integral may be evaluated, and we will express the 
infinitely short propagator as 



P<->(r,i-t / ) 



-iS{t-t') 



Q\\ 
qx -igr 
ig r g± 



(B.12) 



The coefficients are for ao — a\ > 0, given by 
-fe? r -4a + 2ai 4a + 2ai 



Q\\ 



3Traf I ^/a - ai ^/a + 

-.2 q 1 1„2 .),,2 



ai J 



(B.13a) 

1 „2 . 



— fc 3 r 2a 2 , — 3a ai + -^a\ 2a\ + 3a ai + -^a\ 
^ = 3^fl (a - ai) 3 / 2 



} 



(ao + ai)3/ 2 

(B.13b) 



Qr 6™ 2 l(ao- ai ) 3 / 2 (ao + ax) 3 / 2 /- 1 



13c) 



Appendix C: RECIPROCAL EQUATION FOR 
GREEN'S FUNCTION 



In this appendix we derive the reciprocal equation for 
the Green's function. Before doing so we will need some 
results concerning the representation of the Green's func- 
tion. Let us define the following inner product: 



(</)|V) = / d A rdt M(r)(f>{r,t) ■ tf(r,t) 



(C.l) 



We will generally work in the L 2 -space equipped with 
this inner product. Using that the matrix operator M 
is Hermitian, one finds the differential operator T> given 
in equation (jlV.ip to be Hermitian in our inner product 
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space 



(C.2) 



That T> is Hermitian means that the eigenfunctions Fk 
to V 



DF k (r,() = A k F k (r,i), 



(C.3) 



define a complete basis of our inner product space {Fk}. 
A representation of the identity functional given in equa- 
tion (|IV.2j) may therefore be 



£)F£(r,t)F k (r ,to). 



(C.4) 



It can be checked that this is exactly a functional identity 
representation in our inner product space by expanding 
any function on the basis {Fk}. 

To get a formal expression of the Green's function de- 
fined in equation PV.2P we expand the Green's function 
in this basis, and using equation (|C.3[) and (|C4|) we find 



G(r, t\r ,t ) = ^ F k( r ' *)Fk(r , *o). 
k Ak 



(C.5) 



Starting from equation pV.2[) we make the substitution 
t — ► — t, to — ► —t\ and ro — * ri and we write: 



V*G{r, -t\r u -ii) = IS(r, rx)S(t, h). 



(C.6) 



In the next step we take inner product with equa- 
tion pV.2[) and G(r, — t\ri, — t\) from the left with re- 
spect to unprimed coordinates, and equation (|C.6|) and 
G(r, ijro, to) also from the left with respect to unprimed 
coordinates. The resulting two equations are then sub- 
tracted. The term containing u>1 vanish trivially, and 
using rules for differentiating a product, the resulting 
equation may be written as 



2lLU, 



d 

d 3 rdtM(r)— G(r,-i|ri,-*i)-G(r,*|r ,t ) 



d s dt 



M{r)G{Y,-t\v 1 ,-t 1 ) • V x V x M(r) 



G(r,t|r ,t ) - M(r)G(r, i|r , t ) •VxVxM(r) 



G(T,-t\Ti,-h) 



G(ri,ii|r ,t ) - G(r , -io|""i, 

(C.7) 



Using the cut-off property of the Green's function, the 
first term on the left hand side is seen to vanish. Us- 
ing the explicit expression for the Green's function (IC.5j) 
along with Gauss' theorem, one may show that the sec- 
ond term also vanish. The final result is therefore 



G(ri,ti|ro,to) = G(r , -fo|ri,-ii). 



(C.8) 



From Eq. pV.2p . (|C.8|) and using the substitutions t — > 
— t', to — > t, r — > r' and rg — > r we end up with the 



reciprocal equation 
( d o 



c 2 V x V'xM(r'))G(r,t\r',t') 

= I8(r,r')S(t,t'). (C.9) 



In the following we derive the general solution to the 
equation 



d_ 

'dt 



' 2iu},,^- - + c? 



V x V x M(r)Wr,t) = p(r,t), 

(CIO) 



where ^(r, t) is an unknown field, p(r, t) is a source term 
effecting the solution, and M. is some Hermitian matrix 
operator, which may depend on position. We make an 
inner product of equation (|C.10[) with G(r, t\r' ,t') from 
the left and an inner product of equation (|C.9|) with 
ip(r,t) from the right and subtract these two equations. 
In this calculation we are integrating over the time inter- 
val t' € ] to, t + [, where we understand t + — lim e ^o[£ + e]< 
Again we find that terms containing u>1 vanish. Similar 
to above we will use rules for differentiation a product, 
and we eventually end up with 



</>(r,i) 



2mj t 



d 3 r'dt'M{r')G{r, t\r',t') ■ p(r',t') = 

- d r - 
d 3 r'dt'M{r')^- G(r,t\r',t')-ip(r',t') 



dt 



'dt'M(r'){ 
ip(r',t') ■ V x V x M(r')G(r,t\r',t') 

-G(r,t|r',i')- V'x V'x V(r',i')}- 

(C.ll) 



Using the same boundary conditions as was done in the 
calculation leading to the reciprocal equation we conclude 
that the last term in equation (jC.llj) vanish. The right 
hand side of the equation thus reduce to 



fa 



-2iuj L J G?r'M{v') |G(r, t\r',t') ■ ip(r', t') 

2i^J d 3 r'M(r')G(r,t\r',t )-iP(r',to). (C.12) 

Here we have used that the upper time limit vanish due 
to the cut-off in the Green's function. Rearranging terms 
we finally arrive at the general solution to the diffusion 
equation 

</>(r, t) =2i^ J d 3 r' M{r')G(r, t\r', to) ■ tp(r', to) 
+ J£ d 3 r'dt' M(r')&(r,t\r',t') -P{r',t'). (C.13) 
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Appendix D: LORENTZ-LORENZ RELATION 

In the main text we mainly consider lowest order cor- 
rections to the index of refraction. To verify that our 
theory can also correctly reproduce higher order correc- 
tions, we shall in this appendix show how to derive the 
so called Lorentz-Lorenz or Clausius-Mossotti relation for 
the electric permittivity within our theoretical framework 
[37l ]. To lowest order the permittivity is given by Eq. 
(lIlLlll 



elrr^l-V^J]. 



(D.l) 



To calculate the higher order correction it is convenient 
to first Fourier transform the Dyson equation pV.9l) de- 
scribing the light field with respect to time 

(-)/ 



D(-)(r,c) = D{ ) - ) (r^) 



dV P ( -\r,r',uj) ■ m[J]'D ( -)(r',w) 



(D.2) 



This equation is the starting point for the analysis. ( The 
Fourier transformation is here defined as 



(D.3) 



where rj is an infinitely small convergence factor.) 

From Eq. (|IV.24[) we find the Fourier transformed 
propagator P* - ' to read 



*«fk(r') 



(DA) 



The real space representation of this propagator is in gen- 
eral difficult to calculate, however, for a scalar interaction 
the calculation simplify considerably. For uj sa which is 
reasonable in our case, since we are dealing with slowly 
varying operators, the propagator reads 



P (_) (n) = 



d 3 k 



E 



££ 



^,3 gifcL" 

c 2 A-k ki.n 



kV kn 

3i 



[ ( + fc T ,n 



nn 



==- 1 



1 

kj.n 



(D.5) 



where n = r — r', n = |n|, and / is the identity matrix. 
We notice that the propagator gives us the well known 
result for the radiated field of an oscillating dipole. In ad- 
dition we have a term describing a self-interaction. This 
propagator is also discussed in Ref. [35| . In the following 
we shall only be considering the self interaction part of 
the propagator. 

When considering the density correlation function to 
second order (p(ri)p(r 2 )) we have so far used the ideal 
gas approximation in Eq. pV.13[) . where there are no 
correlations between different atoms. In reality we can 
never have two atoms at the same position and this give 
a small correction to (p(ri)p(r 2 )}, which must vanish for 
rx = r 2 (apart from the delta function, which represent 
the single atom contribution). This can formally be de- 
scribed by introducing so called irreducible correlation 
functions h 2 such that 

(p(ri)p(r 2 )) = (p(n)) (p(v 2 )) + h 2 (ri, r 2 ), (D.6) 

where h 2 now takes care of the core-repulsion of the 
atoms (here we exclude the delta function). For ri = r 2 
we thus finds that h 2 (ri,Ti) = —(p(ri)) 2 . 

The above can be used along with the real space repre- 
sentation of the propagator to give the second order cor- 
rection to the permittivity. We will not consider terms 



that vanish when we take quantum mechanical mean. 
The relevant part of the second order term thus gives in 
shorthand notation -/ p(-)(2/3)(V*[J]) 2 E)(-). When 
we introduce this interaction to the differential equation 
(IIII.lip we find the permittivity to second order 



!(r)" 



l_V t [J] + -(V t [J]) 2 . (D.7) 



The calculation can be continued to infinite order [3f 
and the result reads 



itrr^i-v'M-v'Mx; — v 4 [j] 

n=l V 6 



'i + |v*[J]' 



(D. 



This is the Lorenz-Lorenz relation, and we thus see that 
the effect can be included in the theory by dressing the 
spatial mode functions according to the result above. 
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Appendix E: CALCULATIONS OF 
SECOND-ORDER STOKES GENERATOR 



where we have introduced the coefficients 



In this appendix we present detailed calculations of 
the second-order terms of Eq. (|V.31|) . We will denote 
the fourth term of the right hand side of Eq. (|V.31|) as 
, and one finds 



K{$* m3 {v,t)\Sf\r km , r {r,t))) = f-J (fc L /? Cl ) 2 
J f dVdV'^p(r)p(r')e;r(r)*e™,"'(r')4„ i ^n' 



l 

I'ri 



(E.l) 



The seventh term of the right hand side of Eq. (|V.31|) 

=(2) 

plus its complex conjugate we will denote as Sg . To cal- 
culate this term we extend the limits of the time integra- 
tion from minus to plus infinity. This we can do by intro- 
ducing a factor of one half, and approximating the imag- 
inary term i J ^ dt sin(wt) to be zero. This corresponds 
to the usual treatment of such terms in the Markov ap- 
proximation to spontaneous emission when one ignores 
the Lamb shift. We then find the following contribution 
to the Stokes operators 



Kir km3 {v,t)\Sf\r km , 3 ,{v,t))) = (±] (fc L/3ci) 2 



^krn' j' 



+ Qf l n {r)@^'{v')a\ mj a kntv ). 

(E.2) 



In 
I'n 



One notice that the factors of 1 /2 in Eq. (|E.1[) and (|E.2[) 
exactly add up to give one half of the square of the first- 
order term, as is shown in Eq. (IVI.2I) 

The sixth term on the right hand side of Eq. (|V.31[) , 

= (2) 

plus its Hermitian conjugate, we will denote as Sq , and 
we find 



C#"(r) =e,-(r) • {(J(r) x [ e ,,(r) x e ,«(r)]) x e z (r)}. 

(E.4a) 

This term cam be shown to vanish by expanding the spin- 
operator J on the basis defined by the polarization vec- 
tors e x (r), e y (r) and e z (r) and using that the indices 
j, I, I' and I" only run over x and y. 



Finally we will calculate the effect of the fifth term on 
the right hand side of Eq. (|V.31|) , which we will denote 

S D . In this calculation it is important to remember that 
the term will scale as (3 2 p, and reads 



K((r kmj (r,t)\S^\r km>J ,(r,t))) = i^-j (fc L /3) 2 
d 3 r ]T P (r)n m (^T' n '(r)a{ nl a knlll j 

In I 

In 

i(j( r ) • e z( r )) (SjySuc - 5 0X 5iy) {5j> y 5 Vx - 5j' X 5 Vy ) 
+ <%3{r)% l 6 i , l \. (E.5) 



K((r kmj (r,t)\S^\r kmlj ,(r,t))) 



J d 3 r' p(rO{</V)*r(OX V V) 

In 
ql' n' 

if iii 
n I 

®qn , l , ®Q n,/l/, ®knl® km '3 / ^nij^qn'V^ 1 ^'' 1 '' ^ knl 

C#VWV)*? V V)*}, 
(E.3) 



Appendix F: CALCULATION OF 
SECOND-ORDER SPIN-TERMS 



In this section we calculate the second order terms for 
the atomic spin, represented as the third and fourth term 
of the right-hand side of Eq. (|V.42[) . These terms we will 

(2) 

denote J A , and using the previous notation one finds 
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r (2) _-i ( ficik N 1 



=-{—r-\ E / dV E (JW x e *( r )) ( ( ^( r ') ] • e -( r ') J p(OaL'/^-"' 

V ' k J m'm" \ \ J z (r') J J 

E {*r m "(r)*r'" l (r') - n' m (r)*r i "(r')}. (F.l) 

We can examine this term by assuming that the only photon carrying modes of the light are the two modes fk ox and 
iko'y and neglect all other modes. In this case the term reduce to 

j{ 2) = ^p)j 2 J2 J d 3 r' E W x e -( r ))( ( ) •e z (r'))p(r')6L«wIm[*r(r)*r(r')]- 



(n,i)G{(o,a;) J (o' J y)} 



(F.2) 



This term represents an atom at position r' interacting 
with the light field and emitting a photon into mode m, 
which propagates to the position r, where it is absorbed 
by an atom followed by stimulated emission into the clas- 
sical beam. This process is also known as optically in- 
duced dipole-dipole interaction, and indeed the sum over 
all modes m can be used to introduce the dipole propa- 
gator in (|D.5[) . Note, however, that above we have writ- 
ten the term in the paraxial approximation, where we 
ignore the dependence of the polarization vector on the 
mode number. Since the sum over m involves all modes, 
and not just the paraxial modes, an accurate treatment 
requires a more complicated expression involving the po- 
larization vectors along the lines of Appendix [H] (we use 
this more complicated expression in our estimates of the 
size of the effect). 

The last term we will consider is the term describing 
an atom interacting with the light field at two different 
times. This term is represented as the fifth term on the 
right hand side of Eq. (|V.42|1 . and is given on vector 
component form in Eq. (|V. 13|) . We will denote this term 

(2) 

with ' . A short calculation gives 



obtain 



Kfs f mm' 
Tin' 



a kmy a k'm'y aknxak ' n ' x a kmx a k' m' x^ kn y^ k ' n ' V } ' 



At 



(F.4) 



The first order term in Eq. (|V.43[) describe the first order 
effect of rotation of the spin around the e z (r) axis. The 
second order term in (|F.4[) describe the second order term 
of this rotation. From the rotation frequency in the first 
order term cx S3 (assuming '5 to be real), one would thus 
expect this term to scale as /3 2 (s3) 2 which is different 
from the term in (|F.4[) . This difference arises because we 
have separated the term into normal ordered components 
such that the second order term in (IF.4|) only contributes 
when at least two photons are present. When we did 
the normal ordering in the diagram we introduced an 
additional term, which we described by the third term in 
Eq. (IV41) 



7(2) 



1 ( Pcxkt 
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EEE e * ( J - e J 
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Appendix G: CALCULATION OF SPONTANEOUS 
EMISSION 

In this section we calculate the corrections to Eq. 
(|VI.2p and Eq. (|VI.3[) . due to the incoherent interaction. 
To do this we need a result for the infinitely short propa- 
gator. From the definition of the propagator (|IV.24|) and 
the calculation of in it (|V.8|) . we find the relation 



^|f/„(rx)| 2 = ^(r ± ) 



(G.l) 



where we have suppressed the spatial dependence to 
shorten the notation. Doing the sum over and / we 



where g(r±) — kl/(l6ir 2 ) is the zeroth order term of 
the expansion of £>||(r) in (3 given in Eq. (|B.13[) . This 
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= (2) 

result is important when calculating and for relating 
this term with the incoherent interactions, responsible for 
spontaneous emission. When including this term and the 

I 



decay described in Sec. IV B[ the incoherent interaction 
reduce to 



h,out(*±) =■■■- ^ fcL | (r±) / dz'p(z'){c 2 (J 2 (z') - J 2 z (z'))soMr±) + (c 2 J 4 (z') + c 2 [4J 2 (z>) + J y 2 (z')p M „(r ± )}, 

(G.2a) 

s 2 ,out(r±) =■■■- ^ k ^ T±) J dz'p(z')[c 2 3\z') + c 2 [3J 2 (z') + J 2 (z') + J*(z')]}s 2tin (r ± ), (G.2b) 
s 3 ,out(r±)=...- P2Kg 2 (r±) J ^ , p(^){^J 4 (^) + c 1 [J z 2 (z') + ^(^) + ^ 2 (^)]}%n(r±), (G.2c) 



where we have only kept terms that are nonvanishing 
after taking quantum mechanical average of the atomic 
spin. The operators so.m(r_i_) measures the total photon 
flux, and is given as 

so(rx)= \( U *n( r ^)d{ mx a km 'xU m 4r±) 

krnm' 

my a km'y U„Ar±)) (G.3) 



It is important to note that in a discussion of the various 
contributions to decay one should include all terms in 
the perturbative expansion, including the loop diagrams 
(|V.15|) . If these are not included one finds the contri- 
bution from the term in Eq. (|E.5|) to increase the the 
operator S3. 

Similarly we find the effect of spontaneous emission on 
the spin equation to read 



_ ^ 1 

Jx,out(z) = ...- P 2 c 2 1 k L g(r_ L ) {4.«( z )[*o'm( r J + 2*W r -L)] + 7J J y,™( z )sltn( r ±-)} (G.4a) 

fc 

J y ,out(z) =■■■- /3 2 c 2 fc L g(r ± ) ^ { Jy, in {z)[sl in (r ± ) + -s^„(r ± )] + - J x ,i n (z)§l in (r ± )} (G.4b) 

fc 

Jz.out(z) =...- (3 2 c 2 1 k L g(r ± ) ^ J z , ln {z)s k 0ln {r 1 _). (G.4c) 

fc 

I 



The above result is derived from Eq. (|V.9j) by using the 
paraxial approximation and only keeping terms of order 
(3 2 . A minor correction is introduced since we in Eq. 
(|VI.3[) chose a representation that was in fact not normal 
ordered. 



Appendix H: BEYOND PARAXIAL 
APPROXIMATION 

In this section we will go slightly beyond the approxi- 
mation made in Eq. (|V.30[) . and consider the set 

f q (r) = -}=U nq (r)e nj (r). (H.l) 

We will consider the correction this generalization makes 
to the result given i Eq. (|VI.9|) . and therefore de- 



fine spin-components in the local basis given by the set 
(r),e my (r) and 




(H.2) 



for i <G {x,y,z}. These vectors are defined by the fact 
that, e.g., e oy (r) should be transverse and perpendicular 
to the polarization vector arising from the mode function 
^ fc(r)e 0:E (r). e oz is then defined by 
Similarly for the quantum modes m the definition of e mx 
follow from the fact that it should be perpendicular to 
the polarization vector from the mode J7 m fe(r)e mi ,(r). 

With these definitions Eq. (|VI.9|) gives 
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dV p(r')Re[*r(r)]{je ,(r)[e ra (r) • e mx (r)} - J eox (r) [e ox (r) • e„ l2 (r)]} 
PZt =P% + k ^\[^f J d 3 r' p(r')Im[*r(r)]{^e M (r)[e o:c (r) • e mx (r)] - J Sox (r) [e ox (r) • e m2 (r)]}. 



Similarly we find the correction to Eq. (jVI.lOj) to give 



(H.3a) 
(H.3b) 



J out (r) « J in (r) + k^ Cl ^J2 [Mn°(*)}P£> ~ M^TM]-***] {j m (r) x (e oa: (r) x e„ y (r)) }. (H.4) 
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